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SUMMARY 


Modifications were incorporated into the Sikorsky active transmission 
isolation system to provide the capability of utilizing the system as a rotor 
force measuring device. These included; isolator redesign to improve opera- 
tion and minimize friction, installation of pressure transducers in each iso- 
lator and load cells in series with each torque restraint link. Full scale 
vibration tests performed during this study on a CH-53A helicopter airframe 
verified that these modifications do not degrade the systems wide band isola- 
tion characteristics. 

Bench tests performed on each isolator unit indicated that steady and 
transient loads can be measured to within 1 percent of applied load. Individ- 
ual isolator vibratory load measurement accuracy was determined to be i percent 
Load measurement accuracy was found to be independent of variations in all 
basic isolator operating characteristics . 

Full scale system load calibration tests on the CH-53A airframe estab- 
lished the feasibility of simultaneously providing wide band vibration isola- 
tion and accurate measurement of rotor loads. Principal rotor loads (lift, 
propulsive force, and torque) were measured to within 2 percent of applied 
load. This is within the accuracy projected for this system during the pre- 
design study of the Rotor Systems Research Aircraft (Ref. l). The time his- 
tories of transient loads were accurately determined while relative motions 
at the transmission airframe interface were maintained within design tolerances 
of control systems and high speed engine drive shafts. 

In addition to the measurement of principal rotor forces , tests were per- 
formed to establish the feasibility of measuring vibratory rotor forces. It 
was determined that vibratory forces could be measured with an average error 
of 20 percent of applied force at all frequencies where effective isolation was 
provided. The higher error in vibratory load measurement results from the re- 
quirement to adjust measured load cell and isolator loads to include the 
inertia forces resulting from transmission motion. The principal contributor 
to the increased error was found to be the sensitivity of the accelerometers 
utilized in the tests. Evaluation of piezoelectric accelerometers is therefore 
recommended for improved accuracy in future force measurement tests . It is 
anticipated that through the use of the piezoelectric accelerometers, the 
accuracy of vibratory load measurement can be substantially improved. 

As a result of these tests it is concluded that it is feasible - to use a 
rotor force measurement system in determining steady, transient and vibratory 
rotor loads while providing tunable transmission isolation. It is recommended 
that the active transmission isolation system be further developed for' vibra- 
tion control and rotor force measurement in conjunction with slowed rotor opera 
tion. 
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INTRODUCTION 


A basic requirement of the proposed Rotor System Research Aircraft (RSRA) 
is an in flight rotor force measurement system (Ref. l) . In addition, if the 
RSRA is to accomodate certain rotors which operate over a wide range of blade 
passage frequencies, such as the slowed rotor, an active vibration suppression 
system is required to avoid rot or /airframe dynamic resonances. A system of 
this type would also provide tuning for other rotors whose excitation frequen- 
cies fall beyond the operational bands made available through structural tim- 
ing of the airframe. 

Full scale active transmission isolation has been successfully demonstra- 
ted under USAAMRDL Contract DAAJ02-69-C-0101 (Ref. 2). The isolation system 
installed on a CH-53A helicopter consists of three unidirectional hydropneu- 
matic isolators installed vertically at the transmission/airframe interface 
and horizontal rigid links connecting the transmission to the airframe to pro- 
vide torsional restraint. The system exhibits wide band vibration isolation 
at all frequencies above isolator resonances while limiting steady and trans- 
ient relative airframe /transmission deflections within design targets for con- 
trol systems and high speed engine drive shafts . 

Load calibration tests of an individual active isolator performed under 
NASA Contract NAS1-10850 established the feasibility of utilizing the units 
as load measuring devices. 

The objective of this study is to establish, through full scale experi- 
mental ground tests, the feasibility of providing simultaneous wide band active 
transmission isolation and accurate measurement of steady, transient and vibra- 
tory rotor loads . 
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ISOLATOR MODIFICATIONS 


The active isolators as described in Reference 2, were modified and re- 
worked to remove wear resulting from previous operation and to improve their 
accuracy as load measurement devices. Features of the isolator units are illus- 
trated in Figures 1 and 2. Figure 3 shows the isolator modifications. That 
is, the existing sleeves were replaced and the piston head was turned down and 
* replaced with a self-aligning Beryllium-copper head to reduce friction. A 

tight press-fit stainless steel cylindrical sleeve and a ball bushing were in- 
stalled to eliminate the potential of the piston shaft cocking. The piston 
v shaft was built up through chrome plating and ground to eliminate slop. In 

addition, a stainless steel stop and concentric Beryllium- copper self-aligning, 
seal was installed at the base of the piston shaft. Spacers were installed 
between the sleeve and caps to eliminate cocking, caused by torquing the head 
bolts. Hydraulic pressure transducers were installed and the anti-rotation 
device in the servo-linkage was modified to prevent binding during operation. 

ISOLATOR BENCH TESTING 


The modified active isolators underwent load calibration and functional 
bench testing prior to installation on the test vehicle for the full scale sys- 
tem calibration. Steady, vibratory and transient loads were applied to each 
isolator. Effects of variation in isolator precharge, gain, and damping were 
performed to evaluate the effect of said parameters on load calibration. 

Prior to installation in the bench test rig, dry breakout friction was 
evaluated on each isolator. All three isolators exhibited breakout friction 
of under 7 pounds over the limits of piston travel, as is illustrated in Fig- 
ure k. These checks were made during each stage of assembly, to identify 
whether misalignment existed, and also at the conclusion of final assembly. 

The tests were repeated under hydraulic pressure in the bench test rig. No 
change in breakout friction was identified. 

Isolator testing was performed in the bench test facility shown in Fig- 
ure 5- Mounted in series with the instrumented isolator unit was a hydraulic 
ram used for applying steady and vibratory loads and a load cell used to mea- 
sure the applied load. 

Instrumentation, shown in Figure 6, consisted of two absolute pressure 
transducers, a displacement potentiometer for monitoring isolator operation, 
pressure monitoring gages, and three load cells to establish the applied load. 

* Load cells of 5, 10, and 20 thousand pounds were utilized to minimize the 

applied load error in each test range. Oscillograph traces were utilized to 
record direct measurement of applied load, isolator pressures and displacement. 
In addition, two traces were utilized to record the isolator load and error 
which are electronically calculated from direct measurements by a specially 
constructed analog circuit. 
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ISOLATOR BENCH TESTING (continued) 


Basic testing consisted of steady load calibrations including an eval- 
uation of the effects of vibratory load levels which can be anticipated in the 
in-flight environment. The vibratory loading utilized, was a maximum of ±1000 
lbs at sample frequencies of 10, 18.5 and 25 HZ. Transient load tests (15,000 
to 2,000 lbs in 2 seconds) were also conducted. The tests were then expanded 
to evaluate the effect of variations in isolator gain, precharge and damping 
on the load measurement accuracy. The parametric variations tested are listed 
in Table 1. 

The results obtained from these tests are summarized below: 

1. With the isolator at its nominal settings, steady loads were measured 
to within ±1$ of the applied load about a linear bias which can be 
removed through calibration. The presence of vibratory load at fre- 
quencies of 10, 18.5 and 25 HZ, did not affect the accuracy of meas- 
uring steady loads (Figure 7). In addition, all transient test points, 
fall within the established steady error band (Figure 8). 

2. All three active isolators exhibit identical measurement accuracy 
(Figures 7, 9, and 10 ). 

3. Vibratory loads can be measured to within ±h% of applied vibratory 
load, independent of steady load applied or excitation frequency. 
(Figures 11 to 13 illustrate these results). 

4. Isolator load measurement accuracy is independent of variations in 
isolator characteristics. Ten combinations of gain, precharge and 
damping were considered. The resulting steady and vibratory cali- 
bration data is presented in Figures l4 to 21. It is concluded 
from these tests, that parametric variations of isolator parameters 
need not be performed during the full scale system test since no 
effect on load measurement accuracy is detected. 


* 
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FULL SCALE TEST DESCRIPTION 


Test Apparatus and Loading Systems 

The ground test facility used to evaluate the active isolation system 
as a rotor loads measurement system is shown in Figure 22. Principal elements 
, are a "bungee suspension system to simulate free flight, three load application 

systems, for application of steady and transient rotor loads, and a hydraulic 
shalcer mounted at the rotor head to simulate steady-state vibratory rotor loads. 

The structural static CH-53A test vehicle, configured with a MK II tail 
rotor pylon and stabilizer, was weighed and ballasted to a 34,200 pound gross 
weight with a neutral center-of-gravity mass distribution. Concentrated weights 
were used to dynamically simulate the. tail rotor and intermediate gearbox, which 
were not installed on the test vehicle.' The engines were not installed or sim- 
ulated to. simplify the test set-up. 

The bungee suspension system labeled in Figure 23 lifted the entire ve- 
hicle off the deck to simulate free flight by' maintaining all rigid body modes 
below 1 HZ. A rigid ballast (shown in Figure 23) was used to simulate the mass 
of the rotor head and blades and to provide attachments to the suspension sys- 
tem, the external longitudinal loading system, and the shaker. The transmis- 
sion was bolted to a large steel plate, with the isolators suspending the fuse- 
lage through holes cut in the skin of the cabin ceiling (Figure 24). The lower 
isolator attachment was an I-beam frame, weighing approximately 700 pounds, 
bolted to the lower caps of the fuselage frame (Figure 25). The weight of the 
transmission attachment plate, approximately 1300. pounds , was offset by using a 
transmission housing without gears. This permitted the upper body to have a 
mass moment of inertia and center-of-gravity simulating that of the actual vehi- 
cle. System installation was largely dictated by the location of the existing 
CH-53A test vehicle structural members. It must be emphasized that the weight 
and complexity of this installation would not exist if the isolation system were 
integrated into the vehicle's basic design. In-plane and torque restraint 
was provided at the transmission base' by four drag struts located symmetrically 
about station 342 and the fuselage centerline (Figure 24). 

A complete assembly of the test facility is illustrated in Figure 2 6 . 
Simulated steady and transient propulsive force was applied to the rotor head 
with a loading system comprised of a hydraulic cylinder, a 5000 pound BLH load 
cell, an electronic quick release mechanism, and a steel cable attached, in 
series between the rotor head and the aircraft tail structure (Figure 27 ) • 

* Loading was applied to the rotor head by pressurizing the cylinder, thus pro- 

ducing tension in the cable between the rotor head and the tail structure. 
Transient loading was accomplished by activating the quick release mechanism 
which released the applied load instantaneously, applying a step-input load to 
the rotor head. 
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FULL SCALE TEST DESCRIPTION (continued) 


Apparatus a nd Loading Systems (continued) 


Steady main rotor torque loading was provided through a pair of hydrau- 
lic cylinders located so as to apply a couple to the main rotor gearbox. The 
cylinders shown in Figures 2k and 28 attached to the transmission adapter 
plate and to a specially designed steel framework, which in turn was fastened 
to the basic airframe. Applied loads were measured with 20,000 pound BLH load 
cells, located in series with the cylinders. The applied torque was reacted 
by the four in-plane links CFigures 2k and 28). The links, of equal length, 
were located symmetrically about station 3^+2 and the fuselage centerline and 
attached'- between the' transmission plate and the steel loading framework. 

Torque loads were measured with four 20,000 pound BLH load cells, one located 
in series with each link. Design of this apparatus provided for location and 
installation of the torque cylinders , drag links , and steel framework at the 
base of the transmission so as not to alter the capability of the system to 
isolate the aircraft from vertical and in-plane vibratory excitations. Spher- 
ical bearings were incorporated at each attachment point to allow movement of 
the upper body with respect to the basic airframe during isolation tests. 

The vertical transient loading system shown in Figures 25 and 29 consis- 
ted of a hydraulic cylinder attached to the base of the transmission and an 
electronic quick release mechanism fastened to a steel framework which in 
turn bolted to the isolator lower attachment frame. A 10,000 pound BLH load 
cell, in series with the cylinder, measured the applied load. Loading was 
accomplished by pressurizing the cylinder which pulled the transmission toward 
the fuselage. Transients were accomplished by activating the quick release 
mechanism which instantaneously released the applied load, providing a step- 
input load to the upper body. Steady lift variations were also accomplished 
using this system. 

Steady-state vibratory loadings were simulated with a hydraulically 
driven 'unidirectional shaker. The shaker consisted of two counter-rotating 
eccentric masses with adjustable unbalance which produced a unidirectional 
sinusoidal excitation proportional to the square of the rotational speed. 

The shaker was hydraulically powered with a large commercial pump. A manual- 
ly operated flow bypass control valve was used to adjust rotational speed. 

The shaker axis was reoriented to provide vertical, longitudinal, and lateral 
excitations and combinations thereof. 
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Test Instrumentation 


Four measurement systems were used to meet the requirements of these 
tests. The first system, involving a normalizing unit and an x-y-y' plotter 
output display, defined the frequency response of the test vehicle with and 
without the isolators activated. Steady and vibratory load measurement data 
was acquired, utilizing the Dymec Load Measurement System and a 3 6 channel 
CEC oscillograph. A narrow-band FM tape system was used to record the trans- 
ient response real-time histories of the load measurement parameters for trans- 
ient load conditions. The FM tape system was also used as a back-up for mea- 
suring vibratory load conditions. The major instrumentation components and 
master test control panel are shown in Figure 30. Active isolators were in- 
strumented in the same manner as in the load calibration bench tests. Accel- 
erometer locations and their orientation with respect to the test vehicle are 
illustrated by arrows in the structural arrangement shown in Figure 31. Two 
sets of accelerometers were utilized, each set located at the individual upper 
and lower center of gravity of the simulated transmission and rotor head. 

This instrumentation configuration was selected so as to preclude any errors re- 
sulting from relative motion due to flexibility in the shaft. The response 
of the transmission and simulated rotor head, combined with the measured vi- 
bratory load in the isolators and torque restraint links, provide the necessary 
measurement of vibratory forces. 
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VIBRATION ISOLATION SUBSTANTIATION 


The ability of the active isolation system to provide wide-band vibration 
isolation has been demonstrated by the tests described in Reference 2. The 
purpose of these tests was to verify that modifications required for load mea- 
surement did not alter these characteristics. 

The steady-state isolator evaluation was made by comparing frequency 1 

response sweeps from the isolated and hard-mounted test vehicle configurations. 

Sweeps in all three principal directions were conducted from 150 to 1800 CPM 

(2.5 HZ to 30 HZ) with a main rotor unidirectional excitation of 850 pounds > 

achieved at the CH-53A 6P frequency (1110 cpm). 

Due to modifications in the individual units and the overall test vehicle 
configuration, including ballast distribution, system operating pressures were 
found to vary from previous test data. In order to maintain the spring rate 
of 10,000 pounds per inch as specified in Reference 2, air precharge pressures 
were modified as required. The precharge pressure calculations for each unit 
are presented in Appendix I. 

Isolation Substantiation 


The result of selected isolated and unisolated response sweeps are pre- 
sented in Figures 32 to U 3 . These sweeps show a substantial reduction of 
fuselage response due to main rotor excitations in all three principal direc- 
tions , thus substantiating that the modifications required for the load cali- 
bration tests, have not altered the systems isolation characteristics. Effec- 
tive isolation in the inplane directions was obtained from an effective 3P 
through the highest frequency tested (555 CPM to 1800 CPM), while the isola- 
tion present in the vertical direction exists over a slightly narrower band. 
The isolation system under discussion, was designed specifically for the CH- 
53A Helicopter, with an articulated rotor system. Experience with five and 
six bladed rotors of this type indicate that the principal vibratory excita- 
tions are in the inplane direction. Therefore, primary consideration in the 
design (Reference 2) was given to optimizing isolation for these excitations. 
The inplane isolator modes at 300 CPM were located as low as possible to ob- 
tain maximum 6P isolation, while minimizing IP amplification. In the origi- 
nal design of this installation the vertical isolator mode at 650 CPM, was 
constrained by the location of the inplane modes and the geometry of the sys- 
tem. As a result, vertical isolation was not optimized at lower frequencies. 
Since some degree of isolation existed throughout the desired frequency range, 
no attempt to further optimize isolation beyond that achieved in the tests of 
Reference 2 was made. Design considerations required for a wider band of ver- 
tical isolation are discussed in Appendix II. 
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VIBRATION ISOLATION SUBSTANTIATION (continued) 


Isolation Substantiation (continued) 

No response sweep is presented for the pilot-vertical station due to a 
vertical excitation because of an instrumentation malfunction. However, Fig- 
ure 33, which shows pilot-vertical response due to a combined longitudinal/ 
vertical excitation is presented. Since isolation is shown for both the com- 
bined longitudinal/vertical excitation and pure longitudinal excitation (Fig- 
ure 4l), it is concluded that isolation at the pilot vertical station exists 
due to vertical excitation alone. This conclusion is further substantiated in 
Figures 33 through 35 which show reduced vertical response at the other fuse- 
lage locations due to vertical excitations. 
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ROTOR LOADS MEASUREMENT AND CORRELATION 


The principal objective of these tests was to verify the feasibility of 
simultaneously providing active transmission isolation and measurement of rotor 
loads. System calibration tests were performed during which known steady, 
steady-state vibratory and transient loadings were applied to the simulated 
rotor head, through load application systems installed on the test vehicle, 
or by known vehicle weight. Transducers installed on the rotor head , trans- 
mission housing, isolators, and torque restraint links were monitored and re- 
corded to obtain those parameters necessary to measure the applied loads. All 
recorded data was processed using the Computer program described in Appendix 
III. The resulting measured loads are compared in Tables 2 and 3, with the 
known applied loads to establish the accuracy of the isolation system as a 
rotor loads measurement device. 

Conditions investigated included: 

a) Combinations of steady rotor lift, propulsive force and main rotor 
torque . 

b) Combinations of principal rotor vibratory excitation at frequencies 
of 10, 18.5 and 30 HZ, in combination with steady rotor lift. 

c) Steady rotor lift and vertical transient load. 

d) Steady rotor lift and transient propulsive force. 

e) Simulated flight conditions including vibratory loads. 

Magnitudes of all rotor loads applied are representative of levels which 
can be reasonably anticipated in actual flight. 

Steady lift was provided by the bungee suspension system which raised 
the entire test vehicle off the deck to simulate free flight. . Measurement of 
the rotor lift was accomplished by monitoring the loads present in each iso- 
lator. Excellent correlation between applied and measured steady rotor lift 
was obtained. The error in measuring rotor lift, alone and in the presence 
of superimposed steady torque , propulsive force , and vibratory loadings , 
averaged 2 percent of the applied lift. Linearity of the measurement system 
was verified by incrementally varying the lower body weight. The average 
error over the entire range was 2 percent of applied lift. These results are 
presented In Figure U4. 

Steady main rotor torque was applied in the presence of steady lift and 
with combinations of propulsive force and vibratory loads, utilizing the 
torque application system described previously. In all cases, the measured 
torque correlated to within 2 percent of applied torque. 
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ROTOR LOADS MEASUREMENT AND CORRELATION (continued) 


The incremental torque loadings presented in Figure 1+5 show that the torque 
measurement is also linear. 

Propulsive forces were measured with the load cells in the longitudinal 
transmission torque restraint links, to a 1 percent accuracy when measured 
alone. The results of incremental loading, to UOOO pounds , (Figure 1*6 ) , also 
exhibited linearity. However, when propulsive force was measured during the 
simultaneous application of main rotor torque, the error in measured propul- 
sive force increased substantially. This decrease in accuracy is largely 
attributed to deficiencies in the test set-up which includes; large secondary 
inplane loads introduced by the torque application system and the in-plane 
restraint links being configured so as to provide simultaneous measurement of 
main rotor torque, side force, and propulsive force. This coupling and result- 
ant decrease in accuracy was predicted in the ,RSRA pre-design study (Ref. l). 
These test results therefore substantiate the requirement for a rotor balance 
configuration in which a single longitudinal load sensor and two lateral load 
sensors are employed to decouple the measurement of propulsive force and main 
rotor torque. Such a configuration will provide maximum propulsive force 
measurement accuracy although sacrificing side force measurement accuracy. 

The 2 % accuracy obtained in measuring all principal rotor forces is consistent 
with the projections of the RSRA pre-design study (Ref. l) . 

Vibratory loads were applied at frequencies of 10, l8.5 and 30 HZ. Ex- 
amination of oscillograph rolls indicated that the sensitivities utilized did 
not provide adequate resolution of the vibratory signals required for accurate 
phase measurement. Test conditions were repeated, utilizing magnetic tape 
and signals were amplified. Examination of this test data at the frequencies 
of 18.5 and 30 HZ where most effective isolation resulted, indicated that all 
vibratory transducer signals were either at 0° or 180° phase angle with 
respect to the reference force signal. A sample of test data is shown in 
Appendix III - D. However, difficulty in reading phase remained at the 10 HZ 
frequency and is attributed to the closer proximity of the excitation frequen- 
cy to resonances and the resulting modal coupling. This difficulty is identi- 
fied as the source of larger errors in measuring vibratory loads at 10 HZ. 

At the 18.5 and 30 HZ test frequencies, the measurement of applied load 
was strongly dependant upon accelerometer data. The average error in measuring 
total vibratory force at these frequencies is 20 percent (Figure i+7). 

Three vibratory test points at which the error in measuring vibratory load 
varied from 25 - 50 percent are presented in Appendix IV. In these cases, the 
contribution of the accelerometer data to the measurement of the vibratory load 
varied from 65 to 130 percent of the net load. The acceleration levels record- 
ed varied from .06 to .23 g's which is equivalent to a range of 3 to 10 percent 
of the full scale sensitivity of the accelerometers. 
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ROTOR LOADS MEASUREMENT AND CORRELATION (continued) 


It is concluded that the future use of piezoelectric accelerometers are re- 
quired to produce improved accuracy of accelerometer and, therefore, vibra- 
tory load measurement. It is projected that vibratory load measurement accu- 
racy can be substantially improved through the use of these transducers. 

It can also be concluded from the results previously discussed that vib- 
ratory rotor loads can be most accurately measured when the airframe is effec- 
tively isolated. Under this condition, measurement of vibratory loads is pre- 
dominately dependent upon measuring amplitude and phase of transmission accel- 
erations. Future tests should make provisions for providing maximum possible 
resolution of vibratory signals, thus maximizing phase readout and vibratory 
load measurement accuracy. 

Transient load tests were conducted to evaluate the ability of the rotor 
balance system to respond to and measure simulated flight maneuver loads. 
Longitudinal and vertical step input loads of U,000 and 10,000 pounds, respec- 
tively, were applied, and data was recorded on an FM tape system. Data was 
digitized and fed into computer programs (Appendix III-B), where it was pro- 
cessed, and the time histories of applied load, measured load, and relative 
motions between the transmission and airframe were defined. The transient 
data presented in Figures i+8 to 53 represent time histories of load and dis- 
placement over an elapsed time of five seconds , and are computed with- a samp- 
ling rate of 2k0 data points per second. Applied and measured loads were com- 
pared for measurement accuracy and responsiveness of the force measuring sys- 
tem. The graphs presented are plotted individually for ease of comparison. 

The time histories of the applied transient loads were determined by the 
load in the force application device and the response of the simulated trans- 
mission/rotor head. The measured forces were tabulated from transducers in, 
the active isolators and torque restraint links . 

As shown, the active isolation system was able to accurately measure the 
time history of the applied force even though the rate of application far ex- 
ceeds that which would be anticipated in flight. Variance between applied and 
measured load is again attributed to accelerometer sensitivities, which in 
this case plays the predominant role in defining the time history of applied 
load. Of particular significance to in-flight operation is the fact that even 
under the application of these severe transient loads, relative angular and 
vertical motion between the transmission and airframe are well within ;the 
design tolerances of control systems and high speed engine drive shafts. 



CONCLUSIONS 


The testing performed under this program has successfully demonstrated 
the feasibility of providing simultaneous active transmission vibration iso- 
lation and measurement of steady, vibratory and transient rotor loads. As 
a result of these tests, the following conclusions are reached: 

BENCH TEST 

(1) Active isolator units measure steady and transient loads to better 
than 1 percent of applied load about a linear bias which can be 
removed through calibration. 

(2) The presence of vibratory loads do not affect steady load measure- 
ment accuracy. 

(3) Vibratory loads can be measured to within 4 percent of applied 
vibratory load, independent of excitation frequency. 

(4) Load measurement accuracy is independent of variations in basic 
isolator operating characteristics which include servo gain, air 
precharge pressures, and hydraulic damping. 

SYSTEM TEST 

(1) Isolation system modifications required for rotor loads measure- 
ment do not affect isolation characteristics. 

(2) The configuration tested produced significant vibration reduc- 
tions over the entire frequency range, above isolator resonances. 

(3) The average error in measuring steady rotor lift and main rotor 
torque is 2 percent of applied load, independent of the flight 
condition simulated and the presence of vibratory loads . 

(4) The average error. in measuring propulsive force is 1 percent of 
applied load, independent of the presence of vibratory loading. 

(5) The configuration tested which was not originally designed for 
rotor force measurement utilizes transducers which simultaneously 
measure both main rotor torque and propulsive force. This coup- 
ling produces a substantial degradation in propulsive force mea- 
surement error when main rotor torque is applied simultaneously. 

This substantiates the results of pre-design studies (.Ref. l) which 
identified this problem and resulted in a proposed RSRA configuration 
which would decouple propulsive force and main rotor torque mea-' 
surement. 
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CONCLUSIONS (continued) 


At frequencies where effective isolation to all excitations was 
obtained (above 3/rev on the CH-53A) , the average error in meas- 
uring vibratory rotor loads was established as 20 percent of 
applied loads . This error is attributed predominantly to the 
accuracy of accelerometers utilized. It is anticipated that the 
use of piezoelectric accelerometers would substantially improve 
this accuracy. 

The active transmission isolation system is able to accurately 
measure the time histories of transient loads while maintaining 
relative displacements between the transmission and airframe 
within design tolerances of control systems and high speed engine 
drive shafts. 



RECOMMENDATIONS 


It is recommended that any Rotor Force measurement system for the RSRA 
contain the following features : 


(1) A single longitudinal load sensor to decouple the measurement 
of propulsive force and main rotor torque. 

(2) Piezoelectric accelerometers to provide greater sensitivity and 
accuracy in determining transmission and rotor load accelerations 
required for transient and vibratory load measurement. 

(3) Provisions in data acquisition and processing systems to provide 
maximum possible resolution of vibratory signals required for 
accurate phase definition. 

It is also recommended that the active isolator system be developed for 
vibration control and rotor loads measurement on the RSRA for slowed rotor 
operation. The recommended steps for further development are: 


(1) Design analysis of overall flight vehicle/isolation system dyna- 
mic characteristics, including: 

(a) Rotor stability, stress levels, control loads, and ground/ 
air resonance. 

(b) Control system/isolator coupling and stability augmenta- 
tion system. 

(2) Detail Design including: 

(a) Variable isolator operating characteristics. 

(b) Rotor Balance instrumentation. 

(c) Fail safe and reliability considerations. 

(d) Hydraulic supply and temperature control. 

(e) Engine and control systems interface. 

(3) Fabrication and testing on the RSRA including: 


(a) Structural proof tests. 

(b) Tiedown tests to substantiate compatibility of dynamic com- 
ponents . 

(c) Flight tests. 
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Figure 1. - Assembled Active Isolator. 
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Figure 5. - Bench Test Facility. 
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Figure 6. - Isolator Instrumentation 
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Figure 7- - Steady Load Calibration at Basic Isolator Settings - Isolator No. 
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Figure 10. - Steady Load Calibration at Basic Isolator Settings - Isolator No. 






Isolator vibratory, lbs 



0 400 800 1200 1600 2000 2400 

Load cell vibratory, lbs 


Figure 11. - Vibratory Load Calibration at Basic Isolator Settings - 
Isolatox No. 1. 
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Figure 13. - Vibratory Load Calibration at Basic Isolator Settings - 
Isolator No. 3. 
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Figure 14. - Steady Load Calibration With Increased Precharge 
psia.) and Gain Variation - Isolator No. 1. 
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Calibration With Precharge and Damping Variations - Isolator No. 
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Figure 18. - Steady Load Calibration With Precharge Variations - Isolator No. 
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Figure 21. - Vibratory Load Calibration With Precharge Variations - 
Isolator No. 3. 
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Figure 22. - Active Isolator /Rotor Force Measurement Test Facility. 
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Figure 23. - Unidirectional Hydraulic Shaker Installation 







Figure 2U. - Active Isolation Rotor Force Measurement System Installation. 










Figure 26. - Test Facility Illustration 





Figure 27. - Steady and Transient Propulsive Load Application System. 
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Figure 29. - Vertical Transient Loading System. 
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Figure 30. - Master Controls and Instrumentation for Steady and Vibratory Loads. 
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Figure 31. - Accelerometer Locations. 
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Figure 32. - Main Rotor Head-Vertical Response due to Vertical Excitation. 
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Figure 33. - Pilot Vertical Response due to Vertical-Longitudinal Excitation 
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Figure 35. - Aft Cabin Vertical Response Due to Vertical Excitation 
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Figure 36. - Main Rotor Head-Lateral Response Due to Lateral Excitation. 
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Figure 37. - Pilot Head-Lateral Response Due to Lateral Excitation 
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Lateral Response Due to Lateral Excitation. 
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Figure 1*1. - Pilot Head-Longitudinal Response Due to Longitudinal Excitat 
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Figure 1+3. - Aft Cabin Head-Longitudinal Response Due to Longitudinal Excitation 
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Figure 49. - Angular Deflection Between Airframe and Transmission Interface 
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TABLE 1. BENCH TEST ISOLATOR PARAMETERS 


ISOLATOR GAIN 
50 

100 (NOMINAL) 

200 

ISOLATOR LOW SIDE PRECHARGE PRESSURE (PSl) 
150 

300 (NOMINAL) 

500 

ISOLATOR HIGH SIDE PRECHARGE PRESSURE (PSl) 
550 

1100 (NOMINAL) 

1820 

FLUID DAMPING ORIFICE DIAMETER (INCHES) 

.201 

.221 (NOMINAL) 




TABLE II 


TEST SUMMARY STEADY APPLIED VERSUS 
MEASURED LOADS 


Steady lift 

Steady Longitudinal 

Steady torque 

Applied 

Measured 

Applied 

Measured 

Applied 

Measured 







lbs 

lbs 

lbs 

lbs 

ft-lbs 

ft-lbs 







34200 

35070 

— 

— 

— 

— 

24500 

25500 

— 

— 

— 

— 

29200 

29750 

— 

— 

— 

— 

31250 

31750 

— 

— 

— 

— 

34000 

34000 

— 

— 

— 

— 

34200 

34620 

4130 

41 70 

— 

— 



34800 

4180 

4220 

— 

— 



34620 

1013 

1147 

— 

— 



33880 

1987 

2037 

— 

— 



32990 

2894 

2994 

— 

— 



34110 

4129 

4220 

— 

— 



34370 

— 

— 

78800 

81000 



34580 

— 

— 

80600 

83000 



34510 

— 

— 

25500 

25900 



33810 

— 

— 

49000 

48700 

1 

1 

34310 

— 

— 

79200 

79100 

34200 

34360 

5435 

3260 

78300 

80100 



TABLE III - TEST SUMMARY OF VIBRATORY APPLIED VERSUS MEASURED LOADS 


Vibratory loads 

Longitudinal 

Measured 

© 

■ 

© 

a 

n 

+ 1 


O 

O 

7 

o 

5 

270/0° 

20/0° 

o 

o 

o 

CN 

1270/0° | 

530/0° 

o 

00 

CO 

T 

o 

CN 

* 

o 

o 

o 

50/0° 

o 

in 

in 

o 

o 

CO 

1310/0° 

o 

o 

o 

in 

580/98° 

0 

o 

o 

0) 

30/0° 

220/131° 

650/0° | 

60/0° 

Applied 

® 

i/i 

<0 

r 

a 

'q 

+1 


O 

o 

o 

o 

o 

o 

o 

o 

o 

720/0° 

o 

o 

o 

o 

o 

500/0° 

0 

o 

o 

o 

o 

o 

0 

o 

o 

510/0° 

990/0° 

495/0° 

0 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

0/0° 

Lateral 

Measured 

® 

i/i 

a 

£ 

a 

£> 
+ 1 


140/54° 

60/0° 

o 

o 

o 

o 

o 

o 

o 

80/0° 

o 

o 

o 

1900/0.3° 

1350/0° 

0 

o 

s 

o 

CO 

in 

120/19° 

140/0° 

50/0° 

890/1.4° 

1460/0° 

530/0° 

910/42° 

1310/0° 

340/0° 

Applied 

— Ibs/phase 


o 

o 

o 

O 

O 

O 

0 

o 

o 

o 

o 

o 

*0 

o 

o 

o 

\ 

o 

750/0° 

o 

o 

o 

o 

o 

500/0° 

o 

o 

o 

o 

o 

o 

o 

o 

o 

510/0° 

990/0° 

495/0° 

360/0° 

990/0° 

350/0° 

Vertical 

Measured 

— Ibs/phase 


o 

CN 

1 

o 

in 

CN 

935/0° 

730/0° 

280/0° 

o 

o 

\ 

in 

fs 

o 

o 

o 

o 

in 

7 

o 

in 

CO 

o 

o 

o 

0 

o 

o 

o 

CO 

1 

o 

o 

o 

o 

o 

350/0° 

910/108° 

560/0° 

o 

o 

o 

o 

CN 

1420/123° 

1020/0° 

360/0° 

Applied 

—Ibs/phase 


720/0° 

o 

o 

o 

o 

o 

500/0° 

o 

o 

s. 

o 

0 

1 

o 

o 

o 

o 

o 

o 

0 

o 

o 

0 

o 

o 

510/0° 

990/0° 

495/0° 

0 

1 

in 

990/0° 

O 

O 

N. 

in 

0) 

* 

360/0° 

990/0° 

350/0° 

Vibratory 

frequency 



N 

X 


o 

to 

CO 

o 

co 

o 

in 

CO 

o 

CO 

o 

in 

co' 

o 

CO 

O 

in 

co' 

o 

CO 

o 

in 

00 

o 

n 

o 

in 

00* 

o 

00 

Steady torque 

Measured 


M 

£ 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

39300 

41100 

41400 

Applied 


ft-lbs 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

38000 

40200 

40200 

Steady Longitudinal 

Measured 


£ 
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1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1609 

1650 

1684 

Applied 


i/> 

£ 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

CN 

3 

in 

O) 

ID 

in 

00 

«-> 

> 

■o 

<0 

® 

w 

Measured 


£ 


33700 

35120 

33780 

34290 

34280 

34240 

35080 

34980 

34930 

35900 

37200 

36100 

34900 

34900 

34980 

34100 

34120 

34000 

Applied 


£ 


[ 34200 

















34200 
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APPENDIX I 


Active Isolator Precharge Calculations 


In establishing the system precharge pressures, it was determined that the 
individual isolator spring rates would be set to the 10,000 pound per inch value 
utilized during the test of Reference 2. 

The equations for both the high and low pressure side spring rates as a 
function of system operating conditions are developed in Reference 2. It is seen 
that the spring rate at high frequencies is proportional to the operating press- 
ure, piston cross-sectional area, and inversely proportional to the operating 
volume between the diaphram and flow restrictor. This equation is given as: 

2 

Ka,b = Pa,b A a,b 

Yea' ,b ' 

The piston area has not changed, and thus spring rate is proportional to the 
ratio of Pa,b to Vca',b' solely. A sample calculation of the required precharge 
pressures for isolator #1 is presented. Calculations for the remaining units are 
identical in nature. From the data for the original configuration (Table I-l), 
it is seen that the high pressure side spring rate is proportional to the follow- 
ing ratio: 

P b 1100 = 68.8 


The high side operating pressure for isolator #1 is measured at 1300 psi. 
Therefore, the new high pressure side operating volume Vcb' required to maintain 
the same spring rate is : 

Vcb' = 1300 = 18.9 

68.8 


The total available air volume is given as : 

Vtb + Vcb = 14.2 + 26.6 = Uo.8 

o 

The total high pressure side air volume under operating conditions for 
isolator #1 is: 


Vtb + Vcb' = lU.2 + 18.9 = 33.1 

The air pressure under operating conditions is equal to the hydraulic 
fluid operating pressure P^. Therefore, assuming adiabatic expansion, the re- 
quired precharge pressures can be determined from the following relationship: 

P H (Vtb + Vcb Q ) = Pb (Vtb + Vcb') 

Therefore, 

Pu = (1300) (33.1) = 1050 psi 

TUOl 


8 0 



In a similar manner for the low pressure side of isolator if 1: 


Vca' = (420) (6) = 8.4 

300 

Vta + Vca = 7.1 + 26.6 = 33.7 

o 

Vta + Vca' = 7.1 + 8.4 = 15.5 

Therefore, 

P T = (420) (15. 5) = 193 psi 

L (33 7f) 

A summary of all measured and calculated parameters for the three units and the 
basic configuration is presented in Table 1-1. 


TABLE 1-1 

ACTIVE ISOLATOR MEASURED 

AND CALCULATED 

OPERATING 

PARAMETERS 


Original Configuration 

Isolator 1 

Isolator 2 

Isolator 3 

Vca 

26.6 

26.6 

26.6 

26.6 

° 

Vcb 

26.6 

26.6 

26.6 

26.6 

° 

Vta 

7.1 

7.1 

7.1 

7.1 

Vtb 

14.2 

l4 .2 

14.2 

14.2 

Pa 

300 

420 

4oo 

250 

Pb 

1100 

1300 

1150 

900 

Vca’ 

6 

8.4 

8 

5 

Vcb* 

l6 

18.9 

16.7 

13.1 

P H 

810 

1050 

870 

600 

P L 

ll6 

193 

179 

90 
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APPENDIX II 


APPLICATIONS OF ACTIVE ISOLATION TO SLOWED ROTORS 


Considering a reduction of 50 percent in main rotor operating speed, the 
resulting N/Rev frequency would be reduced by a factor of 2 to 555 CPM. At 
this frequency, the degree of vertical isolation is not considered acceptable 
(Figures 32 to 35). Shifting of isolator modes to maintain the desired N/Rev 
isolation can be obtained by lowering the existing spring rate in each isolator. 
By lowering the spring rate and maintaining a constant ratio of N/Rev frequency 
to isolator mode frequencies, the amount of isolation at 555 CPM must be com- 
parable to that obtained with the configuration tested at a forcing frequency of 
1110 CPM. The location of inplane and vertical isolator modes, with the lower 
spring rate and reduced rotor speed are compared with the original spectrum 
obtained at normal operating conditions in Figure (il-l). No amplification at 
1/Rev would result since this frequency is also reduced by 50 percent. 

In order to achieve the modes at the indicated frequencies (Figure II-l), 
the isolator spring rates would have to be reduced from the 10,000 pounds per 
inch tested, to 2500 pounds per inch. Assuming system operating pressures to 
remain constant the required air volumes and precharge pressures required to 
produce this spring rate are calculated as follows : 

From the data of Appendix I- K „ 68.8 =17.2 

a ~T~ 

resulting in a new high pressure side operating volume: 

Vcb ' = 1300 = 75-6 

17.2 

The total available air volume remains at ^0.8 in , however, the new required 
air volume under operating conditions is: 

Vtb + Vcb' = lk. 2 + 75-6 = 89.8 

2 

Since the required air volume exceeds the existing volume , by k9 in , an 
additional supply is required. 

Assuming additional air volume capabilities, including those existing, totaled 
100 in3 , the high air precharge pressure could be calculated from the relation- 
ship derived in Appendix I or: 

P H = (1300) (89-8) = 1168 psi 

100 

In a similar manner for the low pressure side: 

K ^ 50 = 12.5 and 

a r 

Vca* = 420 = 33.6 

12.5 


8k 


Vta + Vca' = 7.1 + 33.6 = 40.7 

Since the total existing available air volume is 33-7, the capability for expan- 
sion is considered. Therefore, an air volume capability totaling 50 in^ is used. 
The low air precharge pressure is found to be: 

P T = (420) (40.7) = 336 psi 


This analysis indicates that available air volumes in the prototype units is 
inadequate by approximately 50 cubic inches for the high pressure side and 20 
cubic inches on the low pressure side. 

In order to expand the existing wide band isolation capabilities, air 
volumes which can be utilized as required would be needed in a universal system. 
A conceptual schematic of such a system is presented in Figure II-2.The load 
measurement capability has been determined to be independant of isolator stiff- 
ness . 
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Figure II-l - Effect of Rotor RPM on Required Isolator Spring Rates 
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Figure II-2 - Proposed Universal Isolation System. 


87 








APPENDIX III 


Computer Program Definition 


APPENDIX III-A 


Sikorsky Aircraft has developed computer programs to analyze data produced 
"by isolation and rotor loads testing. This data may he of two kinds: 

(1) That produced in a steady or vibratory test. The data produced from 
actual test cases is keypunched and input to program Y161+A along with 
geometry and mass property information . From this data, the applied and 
measured forces acting on the system are calculated. Samples of 
actual test cases and a complete program listing are presented. A 
listing of this program is presented in Appendix III-B. 

(2) The data resulting from transient testing . The actual data is recorded 
on analog tape, digitized, and used as a series of time points as 
input to program Y16UE, to calculate the forces in time history form. 

The various program in the analysis package are described as follows: 

Y 16 UA ( 1108 ) 


For steady or vibratory conditions: Takes recorded keypunched data, plus 

geometry and mass data in LOADER format , and isolator calibration curves in 
CLOAD format , and performs calculations to obtain the forces acting on the 
system. 

This program handles the steady and vibratory case on its own. For the 
transient case, a sequence of programs must be run in addition to Yl6UA. The 
sequence is as follows: 

Y 16 UB ( 360 ) 


Initializes disk area and stores basic problem description and dimensions. 
Y164C ( 360 ) 


Takes digital tape content in 9-bit words, validates, removes 'sets' insert- 
ed for reduction purposes and stores the result on disk, where one direct access 
record has data for all gages over 72 points . 

yi6Ud ( 360 ) 


Determines calibration of disk data, further validates and summarizes data 
ready for third-pass program, Y16UE. 

yi6Ue ( 360 ) 


Third-pass program. This takes data from disk using subroutine GETPTS, 
together with LOADER and CLOAD data, and generates time histories of the forces 
on the system. These are written on tape. 
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Y16UF (360) 


This program takes the data recorded on tape, and uses time history plotting 
techniques to produce Benson-Lehner plots. 
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APPENDIX III-B 


COMPUTER PROGRAM LISTING 


THIS DECK WAS PROCESSED BY ENG. COMPUTER SERVICES ON JANUARY 30. 1973. 


11 2345678 


M * • 

OIP 

FOR IRBS.IRBS 

• • • ■ 0 • 


U1005 CD COUNT 
00001 

C 



IRBS0010 

00002 

c 



IRBS0020 

00003 

c 

isolator rb system program . 


IRBS0030 

00004 

c 



IRBS0040 

00005 

c 

ALL QUESTIONS - REFER ART MILTON X 1658 


IRBS0050 

00006 

c 

OR G.C. THOMAS X 2368 


IRBS0060 

00007 

c 



IRBS0070 

00008 

c 



IRBS0080 

00009 

c 



IRBS0090 

00010 

c 



IRBS0100 

00011 

c 

SWITCHES 


IRBS0110 

00012 

c 



IRBSO 120 

00013 

c 

ITYPE = 0 STEADY 


IRBS0130 

00014 

c 

= 1 VIBRATORY 


IRBS0140 

00015 

c 

IRIGID = 0 NON-RIGID 


IR8S0150 

00016 

c 

= 1 RIGID 


IRBS0160 

00017 

c 

IEND = 0 STACK RECORDED DATA 


IRBS0170 

00018 

c 

= 1 STACK LOADER DATA 


IRBS0180 

00019 

c 

= -1 STOP 


IRBS0190 

00020 

c 



IRBS0200 

00021 

c 



IRBS0210 

00022 

c 



IRBS0220 

00023 


DIMENSION NCW >LOC (6) » X < 1000 > . Y ( 1000 ) . Z ( 500 ) 


IRBS0230 

00024 


DIMENSION V(200) 


IRBS0240 

00025 


COMMON /IAIBIG/ I A , IB . IG . JSI2E . IPR INT 


IRBS0250 

00026 

c 



IRBS0260 

00027 


DIMENSION 0(17) .FIM (3) .FIT (3) . ERROR (3) 


IRBS0270 

00028 


DIMENSION R (29,3) . IR (29. 3) 


IRBS0280 

00029 


DIMENSION PHV (3) , PLV ( 3) , PHS(3> . PLS(3> » PHP(3> » PLP(3> , 

IRBS0290 

00030 


1 FHV (3) , FLV (3) , XX(3> . YY(3) . FIV(3) » FIVP(3>. 

IRBS0300 

00031 


2 FITS(3), FLCT (4 ) 


IRBS0310 

00032 


DIMENSION FIVPXI3) 


IRBS0320 

00033 

c 



IRBS0330 

00034 


DIMENSION TITLK12) ,TITL2<12) > NAME 1(2) ,NAME2(2) .DATE ( 2 ) 


IRBS0340 

00035 


DIMENSION TITL3( 12) 


IRBS0350 

00036 

c 



IRBS0360 

00037 


REAL IOXUB1 . IOYUB1 , I0XUB2 . I0YUB2 , IOXLR . I0YLB .IOXUB .IOYUB 

9 

IKBS0370 

00038 


1 IXAC .iyac ,mxa .mya .mza .mxam .myam .mzam 

9 

IRBS0380 

00039 


2 mxaa .myaa ,mzaa 


IRBS0390 

00040 

c 



IRBS0400 

00041 


DIMENSION IERR(3) 


IR8S0410 

00042 


DATA I ERR/ ' ERROR 1ERROR2ERROR3 * / 


IRBS0420 

00043 

c 



IRBS0430 

00044 


DIMENSION ITEMS (29) 


IKBS0440 

00045 


data items/'or tapi tapz flcit flczt flc3t flc4T pih 

P1L 

IRBS0450 

00046 


1P2H P2L P3h P3L CEL1 DEL2 DEL3 ACC1 ACC2 ACC3 

ACC4 

IRBS0460 

00047 


2ACC5 ACC6 ACC7 ACC3 ACC9 ACClO ACC11 ACC12 ACC13 '/ 


IRBS0470 

00048 

c 



IRBS0480 

00049 

c 



IRBS0490 

00050 


9b 


■POJfvjM cn -P m f\) »-• *nrnooa)>vocD«< 40 'ui 4 =oii\JM 


EQUIVALENCE 

EQUIVALENCE 

equivalence 

1 < V C 14)(ZUB1 ) 

2 <V( 18)(XU62 ) 

3 (Vt 22>(IOYUB2> 

4 (V( 26) t ZLB ) 

5 (V( 30) ( YI1 ) 

6 t V < 34)rYI3 ) 

7 (V( 381.XLC2 ) 

8 ( V t 42) »YLC3 ) 

9 (V( 46).ZLC4 ) 

A (V< 50) (ALPHA!) 
B (V( 54) ( ALPHA2) 
C (V( 58) (BETA ) 
D (V( 62) (AH ) 
E C V ( 66) (G ) 

EQUIVALENCE 


EQUIVALENCE 


COMMON /RDUSED/ DR 


DIMENSION RD(29) 



3 


4 


5 


6 


7 

8 



. .C 




. . . n . . . . 







(V( 

1) 

* TYPE 

) ( < V 

2) (RIGID ) 





IR8S0500 

00051 

(V( 

7) 

(DEBUG 

) 







IR8S0510 

00052 











IRBS0520 

00053 

< V { 

11) 

'ftUBl 

) ( ( V 

12) 

t XU81 

( IV 

13) 

»YUb1 

) » 

IRBS0530 

00054 

(vc 

15) 

» IOXUB1 ) > (V( 

16) 

( IOYUB1 

( ( V 

17) 

( WUb2 

) » 

IKBS0540 

00055 

<v< 

19) 

»YUB2 

) ( ( V 

20) 

* ZUB2 

( <V( 

21) 

(I0XUB2) ( 

IRBS0550 

00056 

<V( 

23) 

* v«LB 

) * (V( 

24) 

* XLB 

( t V ( 

25) 

»YLB 

) ( 

IRBS0560 

00057 

(V( 

27) 

» IOXLB 

) » ( V 

28) 

- IOYLB 

( < V 

29) 

rXIl 

) * 

IRBS0570 

00058 

(V( 

31) 

* X 12 

) ( ( V 

32) 

(YI2 

» (V( 

33) 

► X 13 

) * 

IRBS0580 

00059 

(V( 

35) 

* XLC1 

) ( <V( 

36) 

* YLC1 

. < V 

37) 

( ZLCl 

) » 

IR8S0590 

00060 

(V( 

39) 

» YLC2 

) ( <v< 

40) 

* ZLC2 

. I V ( 

41) 

(XLC3 

) * 

IRBS0600 

00061 

<v( 

43) 

(ZLC3 

) ( ( V 

44) 

( XLC4 

( I V 

45) 

(YLC4 

) » 

IR8S0610 

00062 

<V( 

47) 

( XTAP1 

) ( C V 

48) 

» YTAP1 

( C V 

49) 

.ztapi 

) » 

IRBS0620 

00063 

(V( 

51) 

( XTAP2 

) ( ( V 

52) 

( YTAP2 

( I V 

53) 

(ZTAP2 

) » 

IRBS0630 

00064 

(V( 

55) 

»XP 

) ( ( V 

56) 

» YP 

( (VI 

57) 

>ZP 

) » 

IRBS0640 

00065 

(V( 

59) 

* XA 

) » ( V 

60) 

»ya 

( (VC 

61) 

(Za 

> ( 

IRBS0650 

00066 

<v( 

63) 

* AL 

) ( C V 

64 ) 

9 PA 

( I V 

65) 

t PR 

) * 

IRBS0660 

00067 

(V( 

67) 

(C V 

) 







IRBS0670 

00068 

(V( 

70) 

*XA1 

) ( (V( 

71) 

* YAl 

. IV 

72) 

» ZAl 

) » 

IRBS0680 

00069 

(vt 

73) 

» X A2 

) ( C V 

74) 

(YA2 

( IV 

75) 

» ZA2 

) » 

IRBS0690 

00070 

(v( 

76) 

» X A3 

) ( C V 

77) 

* Y A3 

( C V 

78) 

(ZA3 

> ( 

IRBS0700 

00071 

(V( 

79) 

» XA4 

) ( ( V 

80) 

t YA4 

( I V ( 

81) 

( ZA4 

) » 

IR8S0710 

00072 

(V( 

92) 

» XA5 

) ( (V< 

83) 

» YA5 

( (VI 

84) 

t ZA5 

) » 

IRBS0720 

00073 

t V ( 

85) 

* XA6 

) ( C V 

86) 

rYA6 

( I V 

87) 

(ZA6 

) » 

IR8S0730 

00074 

< V c 

88) 

» XA7 

) ( ( V 

89) 

»YA7 

( (Vt 

90) 

(ZA7 

>( 

IRBS0740 

00075 

<v< 

91) 

» X A8 

) * ( V( 

92) 

»YA8 

( I V 

93) 

(ZA3 

) » 

IRBS0750 

00076 

(V( 

94) 

» XA9 

) * ( V ( 

95) 

» Y A9 

( (VI 

96) 

» ZA9 

) » 

IRBS0760 

00077 

(V( 

97) 

» XA10 

) (<v< 

98) 

* YAlO 

. (VI 

99) 

( ZAl 0 

)* 

IRBS077D 

00078 

( V < 100 ) 

» XA11 

) * (V(101) 

» YAH 

( (V ( 102 ) 

» ZA 1 1 

>( 

IRBS0780 

00079 

( V ( 103 ) 

» XA12 

) * ( V ( 104 ) 

» Y A12 

( I V( 105) 

(ZA12 

) > 

IRBS0790 

00080 

( V ( 106 ) 

( XA13 

) ( C V ( 107) 

»YA13 

( (VI 108) 

(ZA13 

) » 

IRBS0800 

00081 

( V ( 109 ) 

» SF 

) * ( V ( 110 ) 

* SFREQ 

( (V(lll) 

(XSH 

) * 

IRBS0810 

00082 

( V ( 1 12) 

(YSH 

) ( ( V( 113) 

» ZSH 

( IV 1 1 14 ) 

» THTLAT ) ( 

IRBS0820 

00083 

( V ( 115) 

(THTLON) 







IRBS0830 

00084 

{ V ( 116) (VERVM ) ( ( V ( 1 17 ) (AINPVM) 





00085 











IRBS0840 

00086 











IRBS0850 

00087 

DR 

# 

tapi 

r TAP2 

(FLC1T (FLC2T ( 




IRBS0860 

00088 

FLC3T ( 

FLC4T 

(P1H 

» P1L » P2H » 




IRBS0870 

00089 

P2L 

t 

P3H 

( P3L 

(DELI (DEL2 ( 




IRBS0880 

00090 

DEL3 (ACC1 

( ACC2 

( ACC3 ( ACC4 ( 




I RBS0890 

00091 

ACC5 (ACC6 

( ACC7 

( ACC8 ( ACC9 ( 




IRBS0900 

00092 

ACC10 ( 

ACCll 

(ACC12 

( ACC13 





IRBS0910 

00093 











IR8S0920 

00094 

) (DR) 









IRBS0930 

00095 

drs 

9 

TAPIS 

t T AP2S 

(FLC1TS(FLC2TS( 




IRBS0940 

00096 

FLC3TS( 

FLC4TS 

( PlHS 

(P1LS (P2HS ( 




IRBS0950 

00097 

P2LS ( 

P3HS 

(P3LS 

(DELIS (DEL2S ( 




IRBS0960 

00098 

DEL3S (ACC1S 

( ACC2S 

( ACC3S ( ACC4S ( 




I RBS0970 

00099 

ACC5S ( 

ACC6S 

( ACC7S 

(ACC8S (ACC9S ( 




IRBS0980 

00100 


95 



1 


1 2 3 4 5 

0 0 0 


6 

0 


7 8 

0 0 U1005 CD COUNT 


C 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 

c 


5 

' 6 

7 

8 
9 
A 
B 
C 
D 
E 
F 
G 
H 

EQUIVALENCE <DRS. 


ACC10S.ACC11S.ACC12S.ACC13S. 

DRV .TAP1V .TAP2V .FLC1TV.FLC2TV. 
FLC3TV.FLC4TV . PlHV .P1LV .P2HV , 
P2LV .P3HV »P3LV .DEL1V .DEL2V , 
DEL3V .ACC1V > ACC2V >ACC3V > ACC4V , 
ACC5V »ACC6 V »ACC7V >ACC8V .ACC9V » 
ACC10V.ACC11V.ACC12V.ACC13V. 

DRP * T AP1P .TAP2P .FLC1TP.FLC2TP, 

FLC3TP.FLC4TP.PlHP fPILP .P2HP » 
P2LP >P3HP .P3LP .DELIP .DEL2P . 
DEL3P . ACC1P .ACC2P .ACC3P > ACC4P . 
ACC5P .ACC6P .ACC7P .ACC8P .ACC9P . 
ACC10P.ACC11P»ACC12P»ACC13P 
R(l.l) ) 


EXTERNAL SELLIN 

READ TITLES > USERS > DATE 

READ ( 5> 10 ) TITL1 »NAME1 • NAME2 > DATE 
10 FORMAT (12A6/6A6) 

READ AND PLOT ERROR CURVE DATA 

JSIZE =1000 

CALL CLOAD(NC.LOC.X.YrZ) 

DO 20 1=1.3 

CALL CPLOT(NCCI) .NC.LOC.X.Y.Z) 

20 CONTINUE 

READ LOADER DATA 

DO 30 1=1.200 
30 V ( I ) =0 . 

40 CONTINUE 

READ (5.50)TITL2 
50 FORMAT (12A6) 

ALPHA1 = ALPHA1 * 57.29578 
ALPHA2 = ALPHA2 * 57.29578 
BETA = BETA * 57.29578 
CALL LOADER (V) 

IDEBUG = DEBUG 
IVERVM = VERVM 
INPVM = AINPVM 


CALCULATE MASS DATA 
WUB = WUB1+WUB2 


IRBS0990 

00101 

IRBS1000 

00102 

irbsioio 

00103 

IRBS1020 

00104 

IRBS1030 

00105 

IRBS1040 

00106 

IRBS1050 

00107 

IRBS1060 

00108 

IRBS1070 

00109 

IRBS1080 

00110 

IRBS1090 

00111 

IRBS1100 

00112 

IRBS1110 

00113 

IRBS1120 

00114 

IRBS1130 

00115 

IRBS1140 

00116 

IRBS1150 

00117 

IRBS1160 

00118 

IRBS1170 

00119 

IRBS1180 

00120 

IRBS1200 

00121 

IRBS1210 

00122 

IRBS1220 

00123 

IRBS1230 

00124 

IRBS1240 

00125 

IRBS1250 

00126 

IRBS1260 

00127 

IRBS1270 

00128 

IRBS1280 

00129 

IRBS1290 

00130 

IRBS1300 

00131 

IRBS1310 

00132 

IRBS1320 

00133 

IRBS1330 

00134 

IRBS1340 

00135 

IRBS1350 

00136 

IRBS1360 

00137 

IRBS1370 

00138 

IRBS1380 

00139 

IRBS1390 

00140 

IRBS1400 

00141 

IRBS1410 

00142 

IRBS1190 

00143 

00144 

00145 

IRBS1420 

00146 

IRBS1430 

00147 

IRBS1440 

00148 

IRBS1450 

00149 

IRBS1460 

00150 


96 


112345678 


XUB 

= (WUBl*XUBl+WUB2*XUB2)/WUB 

IRBS1470 

U1005 CD COUNT 
00151 

YUB 

= IWUB1*YUB1+WUB2*YUB2)/WUB 

IRBS1480 

00152 

ZUB 

= <WUB1*ZUB1+WUB2»ZUB2)/WUB 

IRBS1490 

00153 

IOXUB 

= I0XUB1+I0XUB2 + WUBl*t (YUB1-YUB) »*2 +(ZUBl-ZUB)**2)/386. 

IRBS1500 

00154 

1 

+ WUB2*UYUB2-YUe>**2 + (ZUB2-ZUB > **2 ) /386. 

IRBS1510 

00155 

IOYUB 

= IOYUB1+IOYUB2 + WUB1»( ( XUBl-XUB) **2 +<ZUBl-ZUB>*»2)/386. 

IRBS1520 

00156 

1 

+ WUB2*( (XUB2-XUB)*»2 +(ZUB2-ZUB>**2>/386. 

IRBS1530 

00157 



IRBS1540 

00158 

WAC 

= WUB+WLB 

IRBS1550 

00159 

XAC 

= (WUB*XUB+WLB*XLB)/WAC 

IRBS1560 

00160 

YAC 

= (WUB*YUB+WLB»YL8)/WAC 

IRBS1570 

00161 

ZAC 

= (WUB*ZUB+WLB*ZLB)/WAC 

IRBS1580 

00162 

IXAC 

= IOXUB + IOXLB + WUB * ( ( YUB-YAC ) **2 + (ZUB-ZaC>*»2>/386. 

IRBS1590 

00163 

1 

+ WLB *( (XLB-XAC)**2 + <ZLB-ZAC>**2)/386. 

IRBS1600 

00164 

I YAC 

= IOYUB + IOYLB + WUB *( (XUB-XAC)»*2 + (ZUB-ZAC) **2) /386. 

IRBS1610 

00165 

1 

+ WLB *( (XLB-XAC)**2 + (ZLB-ZAC)**2>/386. 

IRBS1620 

00166 



IRBS1630 

00167 

PRINT TITLE INFORMATION 

IRBS1640 

00168 



IRBS1650 

00169 

WRITE 

(6.60)TITL1.NAME1.NAME2>DATE>TITL2 

IRBS1660 

00170 

60 FORMAT ( 1H1 ,29X , 12A6/30X . 72 ( ' -• ) /// 

IRBS1670 

00171 

1 40X, 

'PREPARED BY '.2A6 // 

IRBS1680 

00172 

2 40X , 

'CHECKED BY '<2A6 // 

IRBS1690 

00173 

3 40X. 

■DATE ' t 2A6 /// 

IRBS1700 

00174 

4 30X. 

12A6/30X.72C-') /) 

IRBS1710 

00175 



IRBS1720 

00176 

PRINT : 

INPUT, MASS ANO GEOMETRY 

IRBS1730 

00177 



IRBS1740 

00178 

WRITE 

(6,520) 

IRBS1750 

00179 

WRITE 

(6, 530) WUB 1 >WUB2, WUB, WLB, wAC , XUB1 , XUB2 , XUB,XLB»XAC , YUBl ,YUB2 IrBS1760 

00180 

1»YUB,YLB,YAC»ZUB1>ZUB2>ZUB,ZLB.ZAC.IOXUB1>IOXUB2,IOXUB>IOXLB>IXAC 

, IRBS1770 

00181 

2I0YUB1 > I0YUB2 > IOYUB . IOYLB . I YAC 

IRBS1780 

00182 

WRITE 

(6,540) 

IRBS1790 

00183 

WRITE 

(6,550) (V(I>, 1=29,34) 

IRBS1800 

00184 

WRITE 

(6,560) 

IRBS1810 

00185 

WRITE 

(6,570) ( V ( I > ,1=35,46) 

IRBS1820 

00186 

WRITE 

(6,580) 

IRBS1830 

00187 

WRITE 

(6,590) (V(I) ,1=47,54) 

IRBS1840 

00188 

WRITE 

(6,600) 

IRBS1850 

00189 

DO 70 

J=1 < 13 

IRBS1860 

00190 

J1 = 

3*(J-l)+70 

IRBS1870 

00191 

J2 = 

3*J +69 

IRBS1880 

00192 

WRITE 

(6,620) J, (V( I ) ,I=J1,J2) 

IRBS1890 

00193 

70 CONTINUE 

IRBS1900 

00194 

WRITE 

(6,610) 

IRBS1910 

00195 

WRITE 

(6.630) <V( I ) ,1=55.58) 

IRBS1920 

00196 

WRITE 

(6,640) XSH, YSH , ZSH 

IRBS1930 

00197 

WRITE 

(6,650) (V(I) ,1=59,61) 

IRBS1940 

00198 

WRITE 

(6»660)PA>AH,PB,AL,G>CV 

IRBS1950 

00199 

WRITE 

( 6, 670 1 SF , THTLON . SFREQ, THTLAT 

IRBS1960 

00200 


97 
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c 

convert angles to radians 

IRBS1970 

00201 


ALPHA 1 = ALPHA1/57. 29578 

IRBS1980 

00202 


ALPHA2 = ALPHA2/57. 29578 

IRBS1990 

00203 


BETA = BETA /57. 29578 

IRBS2000 

00204 

c 


IRBS2010 

00205 

c 

calculate geometry 

IRBS2020 

00206 

c 


IRBS2030 

00207 


0< 1) = YI2 -YA 

IRBS2040 

00208 


D< 2) = YA -YI3 

IRBS2050 

00209 


D( 3) = ZA -ZLC1 

IRBS2060 

00210 


Dt 4) = YUB -YA 

IRBS2070 

00211 


D ( 5) = XA -XII 

IRBS20B0 

00212 


D( 6) = XI2 -XA 

IRBS2090 

00213 


D< 7) = XA -XUB 

IRBS2100 

00214 


D ( 8) = YA -YLC1 

IRBS2U0 

00215 


D< 9) = XA -XLC2 

IRBS2120 

00216 


DUO) = YLC3 -YA 

IRBS2130 

00217 


0(111 = XLC4 -XA 

IRBS2140 

00218 


D(12) = XA -XTAP1 

IRBS2150 

00219 


0(13) = YTAPl-YA 

IRBS2160 

00220 


D(14) = XTAP2-XA 

IRBS2170 

00221 


D(15) = YA -YTAP2 

IRBS2180 

00222 


D 1 16) = XP -XA 

IRBS2190 

00223 


D(17) = ZP -ZA 

IRBS2200 

00224 

c 

print calculated geometry temporarily 

IRBS2210 

00225 


WRITE (6 . 80>D 

IRBS2220 

00226 


80 FORMAT (///* *** 0 *»*•// (5F15. 3/) ) 

IRBS2230 

00227 

c 


IRBS2240 

00228 

c 

READ CARO INPUT OF RECORDED DATA. 29 CARDS .NUMBERED 1 - 29 . 

IRBS2250 

00229 

c 


IRBS2260 

00230 

c 

3 NUMBERS PER CARD 

IRBS2270 

00231 

c 


IRBS2280 

00232 


90 CONTINUE 

IRBS2290 

00233 


READ (5.1001TITL3.IEND 

IRBS2300 

00234 


100 F0RMATU2A6.6X.I2) 

IRBS2310 

00235 


DO 170 1=1.29 

IRBS2320 

00236 


READ (S.110)Jr(IR(I.K>.K=1.3> 

IRBS2330 

00237 


IF (J.NE.I) WRITE (6.120) J» I 

IRBS2340 

00238 


110 FORMAT (15.3110) 

IRBS2350 

00239 


120 FORMAT (> *** ERROR IN INDEX OF RECORDED DATA .'/ 

IR8S2360 

00240 


1 ' INDEX READ '. 15/ 

IRBS2370 

00241 


2 ' INDEX SHOULD BE 15) 

IRBS2380 

00242 

c 


IRBS2390 

00243 

c 

PERFORM VARIOUS DATA CHANGES 

IRBS2400 

00244 

c 

DO 130 K=1 .3 

IRBS2410 

00245 


IRBS2420 

00246 


130 R(I»K)=IR(I.K) 

IRBS2430 

00247 

c 

DIVIDE DEFLECTION DATA BY 1000. 

IN8S2440 

00248 


IF (I.LT.14.0R.I.GT.16) GO TO 150 

IRBS2450 

00249 


DO 140 K=l»3 

IRBS2460 

00250 


98 


1 


c 

c 


c 

c 

c 

c 


c 

c 

c 


c 


c 

c 

c 


c 


c 

c 

c 

c 


2 3 4 5 6 7 8 

o ()•••• , . 





140 R(I,K> = R<I>K)/1000. 

IRBS2470 

UiUU3 CU COUNT 

00251 

150 CONTINUE 

IRBS2480 

00252 

PHASE . DIVIDE BY 100. 

IRBS2490 

00253 

Rtl.3) =-R(I»3)/100. +360. 

IRBS2500 

00254 

DIVIDE 17-29 BY 100. 

IRBS2510 

00255 

IF (I.LT.17) GO TO 160 

IRBS2520 

00256 

RCI.l) = R(I.l)/ 100. 

IRBS2530 

00257 

R( 1 .2) = R(1.2)/ 100. 

IRBS2540 

00258 

160 CONTINUE 

IRBS2550 

00259 

170 CONTINUE 

IRBS2560 

00260 

IRBS2570 

00261 

PRINT RECORDED DATA 

IRBS2580 

00262 

IRBS2590 

00263 


IRBS2600 

00264 

WRITE (6.2201TITL3 

IRBS2610 

00265 

WRITE (6.180) (ITEMS! J> . J. (R ( J. I ) . 1=1 .3) . J=1 .29) 

IRBS2620 

00266 

180 FORMAT(///50X> ’RECORDED DATA ’/50X. 13 < •-’ ) // 

IRBS2630 

00267 

110X, • ITEM’ .13X, ’CHANNEL’ .11X. ’STEADY* .13X. ’VIBRATORY • , 12X. ’ PHASE’ 

/IRBS2640 

00268 

2 //(10X»A6»13X»I2>9X>3(4X.F12.3>4X) ) ) 

IRBS2650 

00269 


IRBS2660 

00270 

STEADY. PICK OFF CORRECT RECORDED DATA . 

IRBS2670 

00271 

DO 190 1=1*29 

IRBS2680 

00272 

IRBS2690 

00273 

190 RD( I ) = R ( I . 1 > 

IR8S2700 

00274 

CONVERT PHASE ANGLES TO RADIANS 

IRBS2710 

00275 

DO 200 1=1.29 

IRBS2720 

00276 

200 R(I.3)=R(I>3)/ 57.29578 

IRBS2730 

00277 

CALCULATE TRUE ISOLATOR LOADS 

IRBS2740 

00278 

IRBS2750 

00279 

FIM(l) = PlH*AH - P1L*AL 

IRBS2760 

00280 

IRBS2770 

00281 

FIMI2) = P2H*AH - P2L*AL 

IRBS2780 

00282 

FIM(3) = P3H*AH - P3L*AL 

IRBS2790 

00283 

DO 210 1=1.3 

IR8S2800 

00284 

IRBS2B10 

00285 

call cipiierr(i) .nc.loc.x.y.sellin.fimid .errorii) .dmi.dm2.om3. 2) 

IRBS2820 

00286 

FIT < I ) = FIM(I) + ERROR ( I ) 

IR8S2830 

00287 

210 CONTINUE 

IRBS2840 

00288 

calculate applied loads at a 

IRBS2850 

00289 

IR8S2860 

00290 


IRBS2870 

00291 

SBETA=SIN(BETA) 

IRBS2880 

00292 

IRBS2B90 

00293 

CBETA=COS(BETA) 

IRBS2900 

00294 

VA = WUB+WLB-DR»SBETA 

IRBS2910 

00295 

DA = OR*CBETA 

IRBS2920 

00296 

MYA = DR*CBETA*0(17) + DR*SBETA*D ( 16) 

IRBS2930 

00297 

SA =0. 

IRBS2940 

00298 

MXA — 0 • 

IRBS2950 

00299 

MZA = 0 • 

IRBS2960 

00300 


99 
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c 


0 

IRBS2970 

Ul 005 CD COUNT 
00301 

c 


IRBS2980 

00302 


WRITE (6,220>TITL3 

IRBS2990 

00303 


220 FORMAT (1h1,29X,12A6/30X»72( ) ) 

IRBS3000 

00304 

c 


IRBS3010 

00305 

c 

PRINT EXTERNAL loads aND TRUf ISOLATOR loads 

IRBS3020 

00306 


write ( 6.680) Va-SA.DA. MXa. MYa, MZA, (FITQ) ,1=1,3) 

IRBS3030 

00307 

c 


IRBS3040 

00308 

c 

calculate measuring system loads at a , 

IRBS3050 

00309 

c 

APPLIED SYSTEM LOADS AT A , 

IRBS3060 

00310 

c 

CALCULATED LOADS AT A 

IRBS3070 

00311 

c 


IRBS3080 

00312 


ALPHA = ALPHA1 

IRBS3090 

00313 


SALPHA=SIN(ALPHA) 

IRBS3100 

00314 


CALPHA=COS(ALPHA) 

IRBS3110 

00315 

c 


IRBS3120 

00316 


VAM = FIT(1)+FIT ( 2 ) +F I T ( 3 ) 

IRBS3130 

00317 


VAA = VA-WUB 

IRBS3140 

00318 


SAM = FLC4T - FLC2T 

IRBS3150 

00319 


SAA = (TAP1-TAP2)*CALPHA +SA 

IRBS3160 

00320 


Dam = FLC1T - FLC3T 

IRBS3170 

00321 


DAA = DA + (TAP1-TAP2) *SALPHA 

IRBS3180 

00322 


MXAM = FIT (2) *D ( 1 ) - FIT (3) *D (2) +(FLC4T-FLC2T)*D(3) 

IRBS3190 

00323 


MXAA = MXA + D(3)*CALPHA*(TAP1-TAP2) - WUB*D<4) 

IRBS3200 

00324 


MYAM = FIT<1)*D<5) "(FIT(2)+FIT(3) ) »D (6) +(FLC3T-FLC1T) *D ( 3) 

IRBS3210 

00325 


MYAA = MYA - WUB*0<7) +D ( 3) *SALPHA* ( TAP2-TAP1 ) 

IRBS3220 

00326 


MZAM = FLClT*D(8) +FLC2T*D (9) +FLC3T*D ( 10 > +FLC4T*D ( 11 > 

IRBS3230 

00327 


MZAA = MZA - TAP1*CALPHA*D(12)-TAP1*SALPHA*D(13) 

IRBS3240 

00328 


1 - TAP2*CALPHA*D( 14) -TaP2*S ALPHA *D (15) 

IRBS3250 

00329 

c 


IRBS3260 

00330 


TAP12 = TAP1-TAP2 

IRBS3270 

00331 


VAC = VAM+WUB 

IRBS3280 

00332 


SAC = SAM- TAP12*CALPHA 

IRBS3290 

00333 


DAC = DAM- TAP12*SALPHA 

IRBS3300 

00334 


MXAC = MXAM- TAP12«CALPHA*D(3)+WUB*D<4) 

IRBS3310 

00335 


MYAC = MYAM+ TAP12*SALPHA*D ( 3) +WU8*D (7) 

IRBS3320 

00336 


MZAC = MZAM+ TAPI *CALPHA*D<12>+ TAP1*S ALPHA *D (13) 

IRBS3330 

00337 


1 + TAP2 *CALPHA*D ( 14 ) + TAP2*S ALPHA *0 (15) 

IRBS3340 

00338 

c 


IRBS3350 

00339 

c 

PRINT OUT RESULTS 

IRBS3360 

00340 

c 


IRBS3370 

00341 


ITYPE=TYPE 

IRBS3380 

00342 


IF (ITYPE.NE.O) GO TO 230 

IRBS3390 

00343 


WRITE (6,690) VAM. SaM, DAM. MXAM, MYAM, MZAM. VAA, SAA. DAA, MXAA, MYAA, 

MZAAIRBS340O 

00344 


1, VAC, SAC, DAC, MXAC »mYAC, MZAC 

IRBS3410 

00345 


IF (IEND.EQ.O) GO TO 90 

IRBS3420 

00346 


IF (IEND.EQ.l) GO TO 40 

IRBS3430 

00347 


IF(IEND.EQ.-l)STOP 

IRBS3440 

00348 


230 CONTINUE 

IRBS3450 

00349 


IRIGID = RIGID 

IRBS3460 

00350 


100 


112345678 


IF (IRIGID.EQ.l) GO to 240 

IRBS3470 

U1005 CD COUNT 
00351 

C 

IRBS3480 

00352 

C NON-RIGID ACCELERATIONS 

IRBS3490 

00353 

C 

IRBS3500 

00354 

C VERTICAL ACCELERATIONS . UPPER MASS OF UPPER BODY 

IRBS3510 

00355 

ZDSUBl = ACC1S 

IRBS3520 

00356 

ZDVU81 = ACC1V 

IRBS3530 

00357 

ZOPUB1 = ACC1P 

IRBS3540 

00358 

C LOWER MASS OF UPPER BODY 

IRBS3550 

00359 

ZDSUB2 = ACC1S 

IRBS3560 

00360 

ZDVUB2 = ACC1V 

IRBS3570 

00361 

ZDPUB2 = ACC1P 

IRBS3580 

00362 

c lateral accelerations 

IRBS3590 

00363 

YDSUB1 = ACC7S 

IRBS3600 

00364 

YDVUB1 = ACC7V 

IRBS3610 

00365 

YDPU81 = ACC7P 

IRBS3620 

00366 

YDSU62 = ACC13S 

IRBS3630 

00367 

YDVUB2 = ACC13V 

IRBS3640 

00368 

YDPUB2 = ACC13P 

IRBS3650 

00369 

C LONGITUDINAL ACCELERATIONS 

IRBS3660 

00370 

XDSUB1 = ACC6S 

IRBS3670 

00371 

XDVUBl = ACC6V 

IRBS3680 

00372 

XDPUB1 = ACC6P 

IRBS3690 

00373 

XDSUB2 = ACC12S 

IRBS3700 

00374 

XDVUB2 = ACC12V 

IRBS3710 

00375 

XQPUB2 = ACC12P 

IRBS3720 

00376 

XDPU81 =XDPUB1 » 57.29578 

IRBS3730 

00377 

YDPUB1 =YDPUB1 » 57.29578 

IRBS3740 

00378 

ZDPUBl =ZDPUBl * 57.29578 

IRBS3750 

00379 

XDPUB2= XDPUB2 * 57.29578 

IRBS3760 

00380 

Y0PU82= YDPUB2 * 57.29578 

IRBS3770 

00381 

ZDPUB2= ZDPUB2 * 57.29578 

IRBS3780 

00382 

C 

IRBS3790 

00383 

C PITCH ACCELERATIONS 

IRBS3B00 

00384 

C 

IR6S3810 

00385 

240 CONTINUE 

IRBS3820 

00386 

IF (IRIGID.EQ.l) GO TO 250 

IRBS3830 

00387 

XCOM1 = ACC2V*SIN (ACC2P) - ACC3V*SIN<ACC3P> 

IRBS3840 

00388 

YCOM1 = ACC2V*C0S (ACC2P) - ACC3V*COS(ACC3P) 

IRBS3850 

00389 

PAUB1P = ATAN2(XCOM1>YCOM1) 

IRBS3860 

00390 

PAUB1V = SORT (XC0M1**2 +YC0Ml**2) /(XA3 -Xa2)*386. 

IRBS3870 

00391 

PAUB1S = ( ACC2S-ACC3S) / ( XA3-XA2) *386. 

IRBS3880 

00392 

250 CONTINUE 

IRBS3890 

00393 

XC0M2 = ACC8V«SIN(ACC8P) - ACC9V*SIN(ACC9P> 

IRBS3900 

00394 

YC0M2 = ACC8V*C0S(ACC8P) - ACC9V*C0S ( ACC9P) 

IRBS3910 

00395 

PAUB2P = ATAN2(XC0M2*YC0M2) 

IRBS3920 

00396 

PAUB2V = SORT ( XC0M2**2 + YCOM2**2)/ (XA9-XA8) *386. 

IRBS3930 

00397 

PAUB2S = (ACC8S-ACC9S)/tXA9-XA8)*386. 

IRBS3940 

00398 

C 

IRBS3950 

00399 

C ROLL ACCELERATIONS 

IRBS3960 

00400 


101 


o o o on o non 


i 


i 

0 


2 3 

0 .... 0 


4 

0 


5 

0 
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C IRBS3970 00401 

IF (IRIGID.EQ.l) GO TO 260 IRBS3980 00402 

XCOM3 = ACC5V*SIN(ACC5P) - ACC4V*SIN(ACC4P) IRBS3990 00403 

YC0M3 = ACC5V*C0S( ACC5P) - ACC4V*C0S(ACC4P> IRBS4000 00404 

RaUbIP = ATAN2(XC0M3.YC0M3) IRBS4010 00405 

RAUBlV = SORT ( XC0M3**2 + YC0M3**2)/ (Ya5-Ya4)»386. IRBS4020 00406 

RAUB1S =(ACC5S-ACC4S)/(YA5-YA4)*386. IRBS4030 00407 

260 CONTINUE IRBS4040 00408 

XC0M4 = ACC11V*SIn(ACC11P)-ACC10V*SIN(ACC10P> IRBS4050 00409 

YC0M4 = ACCllV*COS( ACC11P) -ACC10V*COS ( ACC10P) IRBS4060 00410 

RAUB2P = ATAN2 (XCOM4, YC0M4) IRBS4070 00411 

RAUB2V = SQRT(XC0M4**2 + YC0M4**2)/ ( YA11-YA10) *386. IRBS4080 00412 

RAUB2S = ( ACC11S-ACC10S) / ( YA11-YA10) *386. IRBS4090 00413 

IF ( IDEBUG.EQ.O) GO TO 270 IRBS4140 00414 

WRITE (6,280 ) XCOM1 , YCOM1 IRBS4150 00415 

WRITE (6.2901XCOM2.YCOM2 IRBS4160 00416 

WRITE (6.300 JXC0M3.YC0M3 IRBS4170 00417 

WRITE (6.310JXCOM4.YCOM4 IRBS4180 00418 

270 CONTINUE IRBS4190 00419 

280 FORMATt* XCOMl= «,F12.4.' YC0M1=« »F12.4) IRBS4200 00420 

290 FORMAT ( • XC0M2= »,F12.4,» YCOM2=* »F12.4) IRBS4210 00421 

300 FORMAT ( ' XCOM3= ».F12.4.< YC0M3=* .F12.4) IRBS4220 00422 

310 FORMATC XC0M4= *,F12.4.' YC0M4=' .F12.4) IRBS4230 00423 

00424 

MODIFICATIONS TO Take into account NON-ALIGNMENT OF ACCELEROMETERS. 00425 

00426 

JUST FOR NON-RIGID CASE 00427 

00428 

00429 

IF ( IRIGID.EQ.l) GO TO 271 00430 

00431 

UPPER PART OF UPPER BODY. PITCH. 00432 

00433 

IF (ZA3-ZA2) 2000.3000,2000 00434 

3000 IF(ZA3-ZUB1)2000, 2001.2000 00435 

C 00436 

2000 CONTINUE 00437 

DELTX3 = XA3 - XUB1 00438 

DELTX2 = XUBl- XA2 00439 

DELTZ3 = ZA3 - ZUB1 00440 

DELTZ2 = ZUBl- ZA2 00441 

ANGLE3 = ATAN2 ( DELTZ3 » DELTX3 ) 00442 

ANGLE2 = ATAN2 (DELTZ2.DELTX2) 00443 

RAD3 = SORT (DELTX3**2 + DELTZ3**2> 00444 

RAD2 = SQRT(DELTX2**2 + DELTZ2**2> 00445 

PAUB1V = 386. * SORT ( XC0M1**2 + YC0M1**2) / (RAD2* C0S(ANGLE2) 00446 

* + RAD3* C0SIANGLE3) ) 00447 

PAUB1S = 386. * ( ACC2S - ACC3S) / ( RAD2*COS ( ANGLE2 ) + 00448 

* RAD3*C0S ( ANGLES ) ) 00449 

C 00450 
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2001 CONTINUE 

IF (ZA6-ZUB1 ) 3001(2002.3001 

3001 CONTINUE 

DELTX6 = XA6 - XUB1 
DELTZ6 = ZA6 - ZUB1 
ANGLE6 = ATAN2(DELTZ6»DELTX6> 

RAD6 = SQRT<0ELTZ6*«2 4 DELTX6**2> 

XDSUB1 = ACC6S - (RAD6*PAUB1S*SIN(ANGUE6) ) /386. 

XUUREU = ACC6V*C0S(ACC6F) - (Rad6*PAUB1V*SIN(ANGLE6) > 

« * COS(PAUBIP) / 386. 

XoUIMU = ACC6V*SIN(ACC6P) - (RAD6*PAUB1V*SIN(ANGLE6) ) 

» * SIN(PAUBIP) / 386. 

XDPUBl = AT AN2 ( XDUIMU. XDUREU) *57 .29578 
XOVUBl = S0 RT(XDUReU»* 2 + XDUlMU**2) 

C 

C ROLL 
C 

2002 CONTINUE 

IF1ZA5-ZA4) 2003.3002.2003 

3002 IF(ZA5-ZUB1)2003. 2004. 2003 

2003 CONTINUE 
C 

DELTY5 = YA5 -YUB1 

DELTY4 = YUB1-YA4 

DELTZ5 = ZA5 -ZUB1 

DELTZ4 = ZUB1-ZA4 

ANGLES = ATAN2 (0ELTZ5.DELTY5) 

ANGLE4 = ATAN2IDELTZ4.DELTY4) 

RAD5 = SORT ( DELTy5**2 + DELTZ5**2> 

RAD4 = SORT ( DELTY4**2 + DELTZ4**2> 

RAUB1V = 386. * SORT (XCOM3**2 + YC0M3**2) / (RAD5*C0S(ANGLE5> 

* + RAD4*C0S ( ANGLE4 ) 

RAUB1S = 386. » <ACC5S - ACC4S) / (RAD5*COS(ANGLE5) 

* + RAD4*C0S(ANGLE4) ) 

C 

2004 CONTINUE 

IF CZA7-ZUB1 ) 3003.2005.3003 

3003 CONTINUE 

DELTY7 = YA7 - YUB1 
DELTZ7 = ZA7 - ZUB1 
ANGLE7 = ATAN2 (DELTZ7 .DELTY7) 

RAD7 = SORT ( CELT Y7**2 + DELTZ7**2) 

YDSUB1 = ACC7S + RAD7*RAUB1S*SIN(ANGLE7) / 386. 

YOUREU = ACC7V*C0S ( ACC7P) + (RadT*RAUB1V* SIN(ANGLE7)> 

* * COS(RAUBIP) / 386. 

YDUIMU = ACC7V*SIN(ACC7P> + tRAD7»RAUBlV*SIN(ANGLE7) ) 

* * SIN(RAUBIP) / 386. 

YQPUBl = ATAN2 ( YDUIMU. YDUREU) *57.29578 

YDVUBl = SQRT(YDUREU**2 + YDUIMU **2) 

2005 CONTINUE 


00451 

00452 

00453 
0 0454 

00455 

00456 

00457 

00458 

00459 

00460 

00461 

00462 

00463 

00464 

00465 

00466 

00467 

00468 

00469 

00470 

00471 

00472 

00473 

00474 

00475 

00476 

00477 

00478 

00479 

00480 

00481 

00482 

00483 

00484 

00485 

00486 

00487 

00488 

00489 

00490 

00491 

00492 

00493 

00494 

00495 

00496 

00497 

00498 

00499 

00500 


103 


o o o u 


112345678 

0 0 a 0 0 0 0 U1005 CO COUNT 


00501 

LOWER PART OF UPPER BODY . PITCH. 00502 

00503 

00504 

IF<ZaS-Za8> 2006.400G.2006 00505 

4000 IF<Za9-ZU81)2006.2007.?006 00506 

C 00507 

2006 CONTINUE 00508 

DELTX9 = XA9 - XUB2 00509 

DELTX8 = XUB2- XA8 00510 

DELTZ9 = ZA9 - ZUE‘2 00511 

DEUTZ8 = ZUb2- ZA8 00512 

ANGLE9 = ATAN2(CELTZ9.nCLTX9> 00513 

ANGLES = ATAN2 (DELT Z8 * 0ELTX8 ) 00514 

RAD9 = SQRT<0ELTX9»*2 + DELTZ9**2> 00515 

RAD8 = SQRT(DELTX6**2 + 0ELTz8**2) 00516 

PAUB2V = 386.* SQRT (XC0M2**2 + YCOM2**2) / (RAD8*C0S(ANGLE8) 00517 

* + RAD9*C0S(ANGLE9) ) 00518 

PAUB2S = 336.* (ACC8S-ACC9S) / ( RAD8*C0S ( ANGLE8 ) 00519 

* + RAD9*C0S ( ANGLE9 ) ) 00520 

C 00521 

2007 CONTINUE 00522 

IF ( ZA12 - ZUB2) 4001-2008. 4001 00523 

4001 CONTINUE 00524 

DELX12 = XA12 - XUB2 00525 

DELZ12 = ZA12 - ZUB2 00526 

ANGL12 = ATAN2CDELZ12.DELX12) 00527 

RAD12 = SORT (DELZ12**2 + DELX12**2) 00528 

XDSUb2 = ACC12S - CRAD12 *PAUB2S * SIN(ANGL12) ) / 386. 00529 

XDUREL = ACC12V * COSIACC12P) - (RA012+PAUB2V * SIN(ANGL12)) 00530 

* * C0S(PAUB2P) / 386. 00531 

XDUIML = ACC12V * SIN(ACC12P) - <RAD12*PAUB2V * SIN(ANGL12>) 00532 

* * SIN(PAUB2P) / 386. 00533 

X0PUB2 = ATAN2(XDUIML.XDUREL)*57. 29578 00534 

X0VU&2 = SQRT(XDUREL**2 + XDUIML **2 ) 00535 

C 00536 

C ROLL 00537 

C 00538 

2008 CONTINUE 00539 

IFIZAll- ZA10) 2009.4002.2009 00540 

4002 IFIZAll- ZUB2) 2009.2010.2009 00541 

2009 CONTINUE 00542 

C 00543 

DELY11 = YAH - YUR2 00544 

DELY10 = YUB2 - YA10 00545 

DELZ11 = ZA11 - ZUb2 00546 

0ELZ10 = ZUB2 - ZA10 00547 

ANGL11 = ATAN2<DELZ11.DELY11> 00548 

ANGL10 = ATAN2 (DELZ10 .DELY10) 00549 

RADII = SQRT (DELY 11**2 + DELZ11**2> 00550 
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RADIO = SQRT(DELY10**2 + DELZ10*»2) 


U1005 CD COUNT 
00551 


RAUB2V = 386.* SORT ( XC0M4**2 + YC0M4 **2> / (RADII *COS(ANGLU> 


00552 


* + RADIO *COS(ANGL10) ) 


00553 


RAUB2S = 386.* (ACC11S - ACC10S) / (RADII* COS(ANGLll) 


00554 


* + RADIO* COS(ANGLIO) ) 


00555 

c 



00556 

2010 

CONTINUE 


00557 


IF (ZA13 - ZUB2 ) 4003.2011.4003 


00558 

4003 

CONTINUE 


00559 

C 



00560 


DELY13 = YA13 - YU02 


00561 


DELZ13 = ZA13 - ZU02 


00562 


ANGL13 = ATAN2(DELZ13.DELY13> 


00563 


RAD13 = SORT (DELY 13**2 + DELZ13**2> 


00564 


YDSUB2 = ACC13S + (RAD13* RAUB2S * SINIANGL13)) /386. 


00565 


YDUREL = ACC13V * C0S(ACC13P) + (RAD13* RAUB2V * SIN(ANGL13)) 


00566 


* * C0SIRAUB2P) / 386. 


00567 


YOU I ML = ACC13V * SINIACC13P) + (RAD13* RAUB2V * SIN(ANGL13)> 


00568 


* * SIN ( RAUB2P ) / 386. 


00569 


YDPUB2 = ATAN2 ( YDUIML. YDUREL) *57 .29578 


00570 


YDVUB2 = SQRT(YDUREL**2 + YDUIML **2> 


00571 

C 



00572 

2011 

CONTINUE 


00573 

C 



00574 

271 

CONTINUE 


00575 


RAUB1P =RAUb1P * 57.29578 

IRBS4100 

00576 


PAUB1P =PAUB1P * 57.29578 

IRBS4110 

00577 


RAUB2P= RAUB2P * 57.29578 

IRBS4120 

00578 


PAUB2P= PAUB 2P * 57.29578 

IRBS4130 

00579 

c 


IRBS4240 

00580 


IF (IRIGID.EQ.O) GO TO 320 

IRBS4250 

00581 

C 


IRBS4260 

00582 

c 

RIGIO ACCELERATIONS 

IRBS4270 

00583 

c 


IRBS4280 

00584 


ZDSUB = ACC1S 

IRBS4290 

00585 


ZDVUB = ACC1V 

IRBS4300 

00586 


ZDPUb = ACC1P 

IRBS4310 

00587 


YDSUU = ACC7S 

IRBS4320 

00588 


YOVUb = ACC7V 

IRBS4330 

00589 


YDPUB = ACC7P 

IRBS4340 

00590 


XUSUB = ACC6S 

IRBS4350 

00591 


XDVUB = ACC6V 

IRBS4360 

00592 


XUPUd = ACC6P 

IRBS4370 

00593 


PAUBS = PAUB2S 

IRBS4380 

00594 


PAUBV = PAUB2V 

IRBS4390 

00595 


PAUBP = PAU82P 

IRBS4400 

00596 


RAUBS = RAUB2S 

IRBS4410 

00597 


RAUBV = RAUB2V 

IRBS4420 

00598 


KAUBP = RAUb2P 

IRBS4430 

00599 


XDPUb = XDPUB * 57.29578 

IRBS4440 

00600 
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YDPUB = YDPUB * 57.29578 

IRBS4450 

00601 

ZDPUB = ZDPUB * 57,29576 

IRBS4460 

00602 

RAUBP = RAUBP * 57.29578 

IRBS4470 

00603 

PAUBP = PAUBP * 57.29578 

IRBS4480 

00604 


IRBS4490 

00605 

320 CONTINUE 

IRBS4500 

00606 


IR8S4510 

00607 

PRINT ACCELERATIONS 

IRBS4520 

00608 


IRBS4530 

00609 

IF ( IDEBUG.EQ-0) GO TO 330 

IRBS4540 

00610 

IF ( 1RIGID.EQ.0) WRITE (6. 700 ) ZDSUB1 > ZDVUB1 >ZDPUB1 > ZDSUB2 » ZDVUB2r ZIR8S4550 

00611 

1DPUB2 * YDSU81 » YDVUB1 , YDPUB1 1 YDSUB2 . YDVUB2 » YOPUB2 > XDSUB1 < XDVUel > X0PUIRBS4560 

00612 

2Bl.XDSUB2»XDVUB2iXDPUB2>PAUBlPrPAUBlV»PAUBlS»PAUB2P»PAUB2VfPAUB2S> IRBS4570 

00613 

3RaUB1P>RAUB1V>RAUB1S.RAUB2P.RAUB2V.RAUB2S 

IRBS4580 

00614 

IF (IRIGIO.EQ.l) WRITE (6 .710 ) ZdSUB- ZDVUB > ZDPUB >YDSUB»Y oVUB . YDPUB 

. IR8S4590 

00615 

lXDSUB»XDVUBrXDPUB.pAUBSrPAUBViPAUBP»RAUBS.RAUBV. RAUBP 

IRBS4600 

00616 

330 CONTINUE 

IRBS4610 

00617 


IRBS4620 

00618 

CALCULATE inertia FORCES 

IRBS4630 

00619 


IRBS4640 

00620 

IF (IRIGID.EO-1) GO TO 340 

IRBS4650 

00621 


IRBS4660 

00622 

NON-RIGID 

IRBS4670 

00623 


IRBS4680 

00624 

UPPER MASS OF UPPER BODY 

IRBS4690 

00625 

VIUB1S = WUBl*ZDSUBl 

IRBS4700 

00626 

V IUB1V = WUB1«ZDVUB1 

IRBS4710 

00627 

VIUB1P = ZDPUB 1 + 180. 

IRBS4720 

00628 

SIUBlS = WUBl*YOSUBl 

IRBS4730 

00629 

SIUB1V = WUBl*YDVUBl 

IRBS4740 

00630 

5IUB1P = YDPUB1 + 180. 

IRBS4750 

00631 

DIUB1S = WUB1*XDSUB1 

IRBS4760 

00632 

DIUBiV = WUBl*XDVUBl 

IRBS4770 

00633 

DIUB1P = XDPUB1 + 180. 

IRBS4780 

00634 

RIUB1S = I0XUB1*RAUB1S 

IRBS4790 

00635 

RIUB1V = I0XUB1+RAU61V 

IR8S4800 

00636 

RIUB1P = RAUB1P + 180. 

IRBS4810 

00637 

PIUB1S = I0YU31*PAUBlS 

IRBS4820 

00638 

PIUB1V = I0YUB1*PAUB1V 

IRBS4830 

00639 

PIUB1P = PAUB1P + 180 . 

IRBS4840 

00640 

LOWER MASS OF UPPER BODY 

IRBS4850 

00641 

VIUB2S = WUB2*ZDSUb2 

IRBS4860 

00642 

VIUB2V = WUB2*ZDVU32 

IRBS4870 

00643 

VIUB2P = ZDPUB2 + 180. 

IRBS4880 

00644 

SIUB2S = WUB2*YDSUB2 

IRBS4890 

00645 

SIUB2V = WUB2*YDVUb2 

IRBS4900 

00646 

SIUB2P = YDPUB2 + 180. 

IRBS4910 

00647 

DIUB2S 1 WUB2*XDSUB2 

IRBS4920 

00648 

DIUB2V = WUB2*XDVUr2 

IRBS4930 

00649 

DIUB2P = XDPUB2 + 180. 

IRBS4940 

00650 
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RIUB2S = I0XUB2* RAUB2S 


IRBS4950 

U1005 CD COUNT 
00651 

RIUB2V = I0XUB2* RAUB2V 


IRBS4960 

00652 

RIUB2P = RAUB2P + 180. 


IRBS4970 

00653 

PIUB2S = I0YU82* PAUB2S 


IRBS4980 

00654 

PIUB2V = I0YUB2* PAUB2V 


IRBS4990 

00655 

PIUB2P = PAUB2P + 180. 


IRBS5000 

00656 

ALLOW VIBRATORY MASS To BE USED » NOT STEADY MASS 



00657 

IF(IVERVM.NE.O) VIUB2V=VERVm * ZDVUB2 



00658 

IFIINPVM . NE . 0 ) SIUB2V=AINPVM* YDVUB2 



00659 

IF ( INPVM • NE * 0 ) 0IUB2V=AINPVM* X0VUB2 



00660 



IRBS5010 

00661 

GO TO 350 


IRBS5020 

00662 



IRBS5030 

00663 

UPPER BOOY AS A RIGIO SYSTEM 


IRBS5040 

00664 



IRBS5050 

00665 

340 CONTINUE 


IRBS5060 

00666 

VIUBS = WUB»ZDSUB 


IRBS5070 

00667 

VIUBV = WUB»ZDVUB 


IRBS5080 

00668 

VIUBP = ZDPUB + 180. 


IRBS5090 

00669 

SIUBS = ,VUb*YDSUB 


IRBS5100 

00670 

SIUBV = wub»ydvub 


IRBS5110 

00671 

SIUBP = YDPUB + 180. 


IRBS5120 

00672 

DIUBS = WUB*XDSUB 


IR8S5130 

00673 

DIUBV = WUB*XDVUB 


IRBS5140 

00674 

DIUBP = XDPUB + 180. 


IR8S5150 

00675 

RIUBS = IOXUB*RAUBS 


IRBS5160 

00676 

RIUBV = I0XUB«RAUbV 


IRBS5170 

00677 

RIUBP = RAUBP + 180. 


IRBS5180 

00678 

PIUBS = IOYUB*PAUbS 


IRBS5190 

00679 

piubv = ioyub*paubv 


IRBS5200 

00680 

PIUBP = PAUBP + 180. 


IRBS5210 

00681 

ALLOW VIBRATORY MASS To BE USED 



00682 

IF(IVERVM.NE.O) VIUBV =VERVm » ZDVUB 



00683 

IFIINPVM . NE . 0 ) SIUBV =AINPVM* YDVUB 



00684 

IFIINPVM .NE.O) DIUBV =AINPVM* XDVUB 



00685 



IRBS5220 

00686 

350 CONTINUE 


IRBS5230 

00687 



IRBS5240 

00688 

set up vectors from recorded data 


IRBS5250 

00689 



IRBS5260 

00690 

PHV(l) = P1HV 


IRBS5270 

00691 

PhV (2) = P2HV 


IRBS5280 

00692 

PHVI3) = P3HV 


IRBS5290 

00693 

PLV(l) = P1LV 


IRBS5300 

00694 

PLVI2) = P2LV 


IKBS5310 

00695 

PLV13) = P3LV 


IRBS5320 

00696 

PHS(l) = P1HS 


IRBS5330 

00697 

PHSI2) = P2HS 


IRBS5340 

00698 

PHSI3) = P3HS 


IRBS5350 

00699 

PLS(l) = P1LS 


IP.BS5360 

00700 
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PLS(2> = P2LS 

IRBS5370 

00701 

PLS(3) = P3lS 

IRBS5380 

00702 

PHP(l) = P1HP 

IRBS5390 

00703 

PHP(2) = P2HP 

IR0S54OO 

00704 

PHP (3) = P3HP 

IRBS5410 

00705 

PLPd) = PlLP 

IRBS5420 

00706 

PLP(2) = P2LP 

IRBS5430 

00707 

PLP C 3) = P3LP 

IRBS5440 

00708 


IRBS5450 

00709 

CALCULATIONS per isolator 

IRBS5460 

00710 


IRBS5470 

00711 

DO 360 1 = 1.3 

IRBS5480 

00712 

FHV C I ) = AH * PHV ( I ) 

IRBS5490 

00713 

FLVtl) = AL * PLV ( I ) 

IRBS5500 

00714 

IF (IDEBUS.EQ.O) 60 TO 370 

IRBS5510 

00715 

WRITE (6.36G)AH.AL.PHV(I) ,PLV(I) >FHV(I) »FLV<I> .PHP(I) tPLP(I) 

IRBS5520 

00716 

360 FORMAT!' »*» ' > 8F12.4) 

IRBS5530 

00717 

370 CONTINUE 

IRBS5540 

00718 

XX(I) = FHV(I)*SIM(PHP(I>) - FLV ( I ) *SIN( PLP ( I) ) 

IRBS5550 

00719 

YY ( I ) = FHV(I)*COS<PHP<I> ) - FLV ( I ) *C0S( PLP ( I > ) 

IRBS5560 

00720 

FIV ( I ) = SORT ( XX ( I ) **2 + YY ( I) **2) 

IRBS5570 

00721 

FIVP(I)= ATAN2(XX( I) rYY(I) ) 

IRBS5580 

00722 

FIVPX(I)=FIVP(I) * 57.29578 

IRBS5590 

00723 


IRBS5600 

00724 

CORRECTED STEADY LOADS PER ISOLATOR 

IRBS5610 

00725 


IRBS5620 

00726 

F IM ( I ) = PHSlI)»AH - PLS ( I ) *AL 

IRBS5630 

00727 

CALL CIP( IERR( I> <NC>LOC.X,Y.SELLIN.FIM( I ) » ERROR ( I > .DM1. DM2. DM3 >2> 

IRBS5640 

00728 

F ITS 1 1 ) = FIM ( I ) + ERROR ( I ) 

IRBS5650 

00729 


IRBS5660 

00730 

MEASURING SYSTEM LOADS 

IRBS5670 

00731 


IRBS5680 

00732 

FIT ( I) = FITS(I) + FIV(I)*C0S<FIVP(I) ) 

IRBS5690 

00733 


IRBS5700 

00734 

360 CONTINUE 

IRBS5710 

00735 

IF (IDEBUG.ES.O) GO TO 400 

IRBS5720 

00736 

WRITE (6.390) (XX(I) .YY(I ) . FIV( I ) .FI VPX ( I ) . 1=1 . 3) 

IRBS5730 

00737 

390 FORMAT (2X.4F12. 3) 

IRBS5740 

00738 

400 CONTINUE 

IRBS5750 

00739 

, 

IRBS5760 

00740 

FLCT ( 1 ) = FLC1TS + FLC1TV * COS(FLCITP) 

IRBS5770 

00741 

FLCT (2)= FLC2TS + FLC2TV * C0S1FLC2TP) 

IRBS5780 

00742 

FLCT (3) = FLC3TS + FLC3TV * C0S(FLC3TP> 

IRBS5790 

00743 

FLCT (4)= FLC4TS + FLC4TV * C0SCFLC4TP) 

IRBS5800 

00744 


IRBS5810 

00745 

CALCULATION OF SHAKER FORCE (MEASURED) 

IRBS5820 

00746 


IR8S5830 

00747 

VIUBP = VIUBP / 57.29578 

IRBS5840 

00748 

DIUBP = DIUBP / 57.29578 

IRBS5850 

00749 

SIUBP = SIUBP / 57.29576 

IRBS5860 

00750 
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VIUB1P 
VIUB2P 
DIUB1P 
0IUB2P 
SIUB1P 
SIUB2P 
C VERTICAL 
RIG1 


2 

0 .. 

VIUBlP 

VIUB2P 

DIUB1P 

DIUB2P 

SIUB1P 

SIUB2P 


3 



57.29578 

57.29578 

57.29578 

57.29578 

57.29578 

57.29578 


VIUBV * COS(VIUBP) 

IF(IRIGID.NE.l) RIG1= VIUBlV»C0S(VIUBlP)+VIUB2V*C0S(VIUB2P> 
RIG2 = VIUBV » SIN(VIUBP) 

IF( IRIGID.NE.I) RlG2= VIUB1V*SIN(VIUB1P)+VIUB2V»SIN(VIUB2P) 

IF (IDEBUG.EQ.O) GO TO 410 

WRITE <6,420>RIG1,RIG2»BETA>DRV,DRP> (FIV ( I ) > 1=1 ,3) . (FIVP( I ) . 1= 
410 CONTINUE 


420 FORMAT (/// 


RIGl 

— 

, F12.4 

/ 

1 


RIG2 

— 

, F12.4 

/ 

2 


beta 

~ 

, F12.4 

/ 

3 


DRV 

~ 

, F12.4 

/ 

4 


DRP 

~ 

, F12.4 

/ 

5 


FIV 


, 3F12.4 

/ 

6 


FIVP 

~ 

,3F12.4 

) 


SBETA = SIN(BETA) 

FS2R = DRV*SBETA*COS(DRP) - RIGl 

1 + FIV(l)*COS(FIVPtl) )+FIV(2)*C0S(FIVP(2) ) 

2 + FIV (3) *C OS (FI VP (3) ) 

FSZI = DRV*S8ETA*SIN(DRP) - RIG2 

1 + FIV(1)*SIN(FIVP(1) >+FIV(2)*SIN(FIVP(2) ) 

2 ♦ FIV(3)*SIN(FIVP(3) ) 

FSZ = SORT (FSZR**2 + FSZI**2> 

PHIFSZ= ATAN2(FSZI,FSZR) 

C LONGITUDINAL 

RIGl = QIUBV * COS(DIUBP) 

IF (IRIGID.NE.I) RIG1=DIUB1V*COS(DIUB1P)+DIUB2V*COS(DIUB2P) 
RIG2 = DIUBV » SlN(OIUBP) 

IF ( IRIGID.NE.I ) RIG2=DIUB1V*SIN(DIUB1P)+DIU82V*SIN(DIUB2P> 
C 

CBETA = COS(BETA) 

SALPHA= SIN(ALPHA) 

C 

FSXR =-DRV*CfeETA«COS(DRP) - RIGl 

1 -TAP1V*SALPHA»C0S(TAP1P) + TAP2V*SALPHA*C0S(TAP2P) 

2 -FLC3TV*C0S(FLC3TP) + FLC1TV*C0S(FLC1TP) 

FSXI =-DRV*CBETA*SIN(DRP> - RIG2 

1 -TAP1V*SALPHA*SIN(TAP1P) + TAP2V*SALPHA*SIN(TAP2P) 

2 -FLC3TV*SIN(FLC3TP) + FLC1TV*SIN(FLC1TP> 

FSX = SORT (FSXR**2 + FSXI**2) 

PHIFSX= ATAN2( FSXI, FSXR) 

C LATERAL 

RIG2 = SIUBV « SlN(SIUBP) 


8 

U1005 CD COUNT 

IRBS5870 

00751 

IRBS5880 

00752 

IRBS5890 

00753 

IRBS5900 

00754 

IRBS5910 

00755 

IRBS5920 

00756 

IRBS5930 

00757 

IRBS5940 

00758 

IRBS5950 

00759 

IRBS5960 

00760 

IRBS5970 

00761 

IRBS5980 

00762 

IRBS5990 

00763 

i IRBS6000 

00764 

IRBS6010 

00765 

IRBS6020 

00766 

IRBS6030 

00767 

IRBS6040 

00768 

IRBS6050 

00769 

IRBS6060 

00770 

IRBS5070 

00771 

IRBS6080 

00772 

IRBS6090 

00773 

IRBS6100 

00774 

IRBS6110 

00775 

IRBS6120 

00776 

IRBS6130 

00777 

IRB&6140 

00778 

IRBS6150 

00779 

IRBS6160 

00780 

IRBS61 70 

00781 

IRBS5180 

00782 

IRBS6190 

00783 

IRBS6200 

00784 

IRBS6210 

00785 

IRBS6220 

00786 

IRBS6230 

00787 

IRBS6240 

00788 

IRBS6250 

00789 

IRBS6260 

00790 

IRBS6270 

00791 

IRBS6280 

00792 

IRBS6290 

00793 

IRBS6300 

00794 

IRBS6310 

00795 

IRBS6320 

00796 

IRBS6330 

00797 

IRBS6340 

00798 

IRBS6350 

00799 

IRBS6360 

00800 
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IF(IRIGID.NE.I) RIg2=SIUB1V*SIN(SIUB1P>+SIUB2V*SIN(SIUB2P> 

IRBS6370 

U1005 CD COUNT 
00801 

RIGl = 

SIUBV * COS(SIUBP) 

IRBS63B0 

00802 

IP ( IRIGID.NE • 1 ) RIG1=SIUbiV*COS(SIUB1P)+SIUB2V*COS(SIUB2P> 

IRBS6390 

00803 

C 


IRBS6400 

00804 

IF ( IDEBUG»eQ*0) GO to 430 

IRBS6410 

00805 

write (6 

f 440) RIGl >RIG2f TAPI V* ALPHA »TAPlPfTAP2V»TAP2p.FLc2TV 

FLC2TPIRBS6420 

00806 

1.FLC4TV. 

FLC4TP 

IRBS6430 

00807 

430 CONTINUE 


IRBS6440 

00808 

440 FORMAT (///' RIGl = '.F12.4/ 

IRBS6450 

00809 

1 

» RIG2 = * *F12 .4/ 

IRBS6460 

00810 

2 

• TAPIVz '.F12.4/ 

IRBS6470 

00811 

3 

• ALPHA= '.F12.4/ 

IRBS6480 

00812 

4 

' TAPIPr 1 » F12 , 4/ 

IRBS6490 

00813 

5 

• TAP2V= '.F12.4/ 

IRBS6500 

00814 

6 

• TAP2P= '.F12.4/ 

IRBS6510 

00815 

7 

» FLC2TV= »rF12.4/ 

IRBS6520 

00816 

8 

• FLC2TP= '.F12.4/ 

IRBS6530 

00817 

9 

• FLC4TV= '.F12.4/ 

IRBS6540 

00818 

A 

• FLC4TP= '.F12.4 ) 

IRBS6550 

00819 

CALPHA= 

COS(ALPHA) 

IRBS6560 

00820 

C 


IRBS6570 

00821 

FSYR = 

-TAPlV*CALPHA*COS ( TAPlP) + TAP2V»CALPHA*C0S(TAP2P) 

IRBS6580 

00822 

1 

-FLC2TV*C0S<FLC2TP) 4 FLC4TV*C0S(FLC4TP) 

IRBS6590 

00823 

2 

-RIGl 

IRBS6600 

00824 

FSYI = 

-TAP1V*CALPHA»SIN( TAPlP )+ TAP2V*C ALPHA*S I N < TAP2P ) 

IRBS6610 

00825 

1 

-FLC2TV*SIn<FLC2TP> + FLC4TV*SIN(FLC4TP) 

IRBS6620 

00826 

2 

-RIG2 

IRBS6630 

00827 

FSY = 

SORT (FSYR**2 + FSYI**2) 

IRBS6640 

00828 

phifsy= 

AT AN2 (FSYI ,FSYR> 

IRBS6650 

00829 

IF < IDEBUG.EQ.O) GO TO 450 

IRBS6660 

00830 

WRITE (6 

. 460 ) FSZR . FSZ I . FSXR . FSXI » FSYR . FSYI 

IRBS6670 

00831 

450 CONTINUE 


IRBS6680 

00832 

460 FORMAT (///' FSZR = F12.4 / 

IRBS6690 

00833 

1 

' FSZI = • r F12.4 / 

IRBS6700 

00834 

2 

' FSXR = ' . F12.4 / 

IRBS6710 

00835 

3 

• FSXI = ' » F12.4 / 

IRBS6720 

00836 

4 

• FSYR = ' . F12.4 / 

IRBS6730 

00837 

5 

• FSYI = F12.4 ) 

IRBS6740 

00838 

VIUBP = 

VIUBP * 57.29578 

IRBS6750 

00839 

DIUBP = 

DIUBP * 57.29578 

IRBS6760 

00840 

SIUBP = 

SIUBP * 57.29578 

IRBS6770 

00841 

VIUB1P = 

VIUB1P * 57.29578 

IRBS6780 

00842 

VIUB2P = 

VIUB2P * 57.29578 

IRBS6790 

00843 

DIUB1P = 

DIUB1P * 57.29578 

IRBS6800 

00844 

DIUB2P = 

DIUB2P * 57.29578 

IRBS6810 

00845 

SIUB1P = 

SIUB1P * 57.29578 

IRBS6820 

00846 

SIUB2P = 

SIUB2P * 57.29578 

IRBS6830 

00847 

THTLON = 

THTLON / 57.29578 

IRBS6840 

00848 

THTLAT = 

THTLAT / 57.29578 

IRBS6850 

00849 

PHIFSZ = 

PHIFSZ * 57.29578 

IRBS6860 

00850 
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1 1 




2 

a 


4 

o 


5 

o 


6 

0 


7 8 

0 0 U1005 CD COUNT 


PHI FSX = PHIFSX 
PHIFSY = PHIFSY 


57.29576 

57.29576 


C CALCULATION OF SHAKER FORCE (APPLIED) 

C 

FSZA = SF*C0S(THTtON)*COS(THTLAT) 

FSYA = SF*SIN(THTLAT) 

FSXA = SF*SIN(THTLON) 

FIFSZAr 0. 

FIFSYA= 0. 

FIFSXA= 0. 

IF (ITYPE.EG.O) GO TO 510 
C 

C VIBRATORY PRINTOUT 

C 

WRITE (6.720) 

IF ( IRIGID.EO. 1) GO TO 470 
C 

C NON-RIGID ACCELERATIONS 

C 

WRITE (6.730) XDSUB1. YOSUB1 .ZDSUBl .RAU81S.PAUB1S. XDVUB1 
1B1.RAUB1V .PAU81V * XDPUB1 » YoPUBl .ZDPUB1 .RAUB1P.PAUB1P 
WRITE (6.740) XDSUB2. YDSUB2 . ZDSUB2 . RAUB2S. PAUB2S . XDVUB2 
1B2.RAUB2V . PAUB2V . XOPUB2 . YDPUB2 . ZDPUB2 > RAUB2P . PAUB2P 
GO TO 4eo 
470 CONTINUE 
C 

C RIGID ACCELERATIONS 

c 

WRITE ( 6 . 750 ) XDSUB . YOSUB . ZOSUB . R AUBS . PAUBS . XDVUB . YD VUB 
1 . PAUBV . XDPUB * YOPUB . ZDPUB . R AUBP . PAUBP 
480 CONTINUE 
C 

WRITE (6.760) 

IF(IVERVM.NE.O.OR.INPVM.NE.O) WRITE <6. 761) 

761 FORMATdOX. • (VIBRATORY MASS USED)' / ) 

IF ( IRIGID.EQ. 1 ) GO TO 490 
C 

C NON-RIGID INERTIA FORCES 
C 

WRITE (6.730)DIUB1S.SIUB1S.VIUB1S.RIUB1S.PIUB1S.DIUB1V 
11V.RIUB1V.PIUB1V.DIUB1P.SIUB1P.VIUB1P.RIUB1P.PIUB1P 
WRITE (6.740)DIUB2S.SIU62S.VIUB2S.RIUB2S.PIUB2S.0IUB2V 
12V.RIUB2V.PIUB2V.DIUB2P.SIUB2P.VIU82P.RIUB2P.PIUB2P 
GO TO 500 
490 CONTINUE 

RIGID INERTIA FORCES 

WRITE <6.750>OIUBS»SIUBS.VIUBS.RIUBS»PIUBS.DIUbV.SIUBV 


IRBS6870 
IRBS6880 
IRBS6890 
IR8S6900 
IRBS6910 
IRBS6920 
IRBS6930 
IRBS6940 
IRBS6950 
IRBS6960 
IRBS6970 
IRBS6980 
IRBS6990 
IRBS7000 
IRBS7010 
IRBS7020 
IRBS7030 
IRBS7040 
IRBS7050 
IRBS7060 
.YDVUB1.ZDVUIRBS7070 
IRBS7080 
. YDVUB2. ZDVUIRBS7090 
IRBS7100 
IRBS7110 
IRBS7120 
IRBS7130 
IRBS7140 
IRBS7150 
. ZDVUB'RAUBVIrbS7160 
IRBS7170 
IRBS7180 
IRBS7190 
IRBS7200 


IRBS7210 
IRBS7220 
IRBS7230 
IRBS7240 
•SIUB1V . VIUBIRBS7250 
IRBS7260 
.SIUB2V.VIUBIRBS7270 
IRBS72B0 
IRBS7290 
IRBS7300 
IRBS7310 
IRBS7320 
IRBS7330 
> VIUBV .RIUBVIRBS7340 


00851 

00852 

00853 

00854 

00855 

00856 

00857 

00858 

00859 

00860 
00861 
00862 

00863 

00864 

00865 

00866 

00867 

00868 

00869 

00870 

00871 

00872 

00873 

00874 

00875 

00876 
00677 

00878 

00879 

00880 
00881 
00882 

00883 

00884 

00885 

00886 

00887 

00888 

00889 

00890 

00891 

00892 

00893 

00894 

00895 

00896 

00897 

00898 

00899 

00900 


111 


1 


1 

0 


2 

0 


4 

0 


5 

0 


6 7 8 
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C 

C 

C 

c 

c 

c 


c 


c 

c 


l»PIUbVrDIU8P»SIUbP » VIU8 P»RIUBPfPIU3P IRBS7350 

500 continue irbS736o 

IRBS7370 

ISOLATOR loads IRBS7380 

IRBS7390 

WRITE (6.770>FITS#fIV.FIVPX IRBS7400 

IRBS7410 

VIBRATORY OUTPUT IRBS7420 

IRBS7430 

WRITE (6»780) IRBS7440 


WRITE (6f790)DAM»SaM> VAM.MXAM'MYAM f MZAM • DAA # SAA r VAA . MXA A > MYAA > MZAAIRBS7450 


1 * DAC t SAC . VAC > MXAC > MYAC > MZAC IRBS7460 

WRITE (6r800> IRBS7470 

WRITE (6f810)FSZA»FIFSZA»FSYA»FIFSYA»FSXAtFIFSXA»FSZ>PHIFSZ,FSY.PHlRBS7480 


1IFSY fFSXfPHIFSX 

510 CONTINUE 

IF UEND.EG.O) GO TO 90 
IF (IEND.E0.1) GO TO 40 
IF(IEND.EO.-1)STOP 


520 FORMAT(///60X> 'MASS DATA ' /60X ’ 9 ( ' -• ) /// 

1' UPPER PAST OF LOWER PART OF 

2L UPPER BODY LOWER BODY TOTAL AIRCRAFT 

3' LOWER BODY LOWER BODY 

530 FORMAT ( 10X » ' W ' »8X r 5 ( 4Xf F12.3 ' 4X> / 

1 IOXf'X • »8Xf5(4X.F12.3>4X> / 

2 lOXr'Y ' »8X.5<4X.F12,3>4X> / 

3 lOX.'Z • .8X,5(4X.F12,3f4X) / 

4 10X. ' IX' »8X.5(4X»F12.3>4X) / 

5 10X»'IY',8X,5(4X.F12.3i4X) / ) 

540 FORMAT (///60X» ' GEOMETRY */60Xf8('-»)// 

1 30X» 'NO. • rl5X» 'X' .20X. 'Y' .19X, 'Z' . 17Xi 'ANGLE' // 

2 iox»' isolators:'/ > 

550 FORMAT ( 30X»'1'»9X» 2(4X.F12.3,4X) / 

1 30X»'2'>9X, 2(4X»F12.3,4X) / 

2 30X»'3«r9X, 2(4X»F12.3,4X) / ) 

560 format:// IOXi 'LOAD CELLS (MEASURING) :’/ ) 

570 FORMAT ( 30X . • 1 ' > 9X . 3<4X ( FI2.3 r 4X ) / 

1 30X. '2' >9X,3(4XeF12.3.4X) / 

2 30X. '3* .9X.314X.F12. 3. 4X) / 

3 30X. '4' ,9X,3(4X.F12.3»4X) / ) 

580 FORMAT:// lOXt'LOAD CELLS (TORQUE APP):'/ ) 

590 FORMAT: 30X, '1' ,9X.4(4X.F12.3.4X> / 

1 30X> *2' #9Xr4(4X»F12.3*4X) / ) 

600 format:// 10X» 'ACCELEROMETERS: • >22X» 'X' »20X> 'Y* .19X» 'Z' >17x» 
1 'ANGLE' /) 

610 FORMAT:// 10X. 'PROPULSIVE force:' /) 

620 FORMAT: 29X.I2»9X»3(4X»F12.3f4X) ) 


IRBS7490 

IRBS7500 

IRBS7510 

IRBS7520 

IRBS7530 

IRBS7540 

IRBS7550 

IRBS7560 

IRBS7570 

TOTAIRBS7580 

•//IRBS7590 

'/1IRBS7600 

IRBS7610 

IRBS7620 

IRBS7630 

IRBS7640 

IRBS7650 

IRBS7660 

IRBS7670 

IRBS7680 

IRBS7690 

IRBS7700 

IRBS7710 

IRBS7720 

IRBS7730 

IRBS7740 

IRBS7750 

IRBS7760 

IRBS7770 

IRBS7780 

IRBS7790 

IRBS7800 

IRBS7810 

IRBS7820 

IRBS7830 

IRBS7840 


00901 

00902 

00903 

00904 

00905 

00906 

00907 

00908 

00909 

00910 

00911 

00912 

00913 

00914 

00915 

00916 

00917 

00918 

00919 

00920 

00921 

00922 

00923 

00924 

00925 

00926 

00927 

00928 

00929 

00930 

00931 

00932 

00933 

00934 

00935 

00936 

00937 

00938 

00939 

00940 

00941 

00942 

00943 

00944 

00945 

00946 

00947 

00948 

00949 

00950 
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0 
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630 FORMAT ( 40X. 4<4X.F12.3.4X1 /) 

640 FORMAT!// 10X. 'SHAKER FORCE:' / 

1 40X> 3(4X.F12. 3.4X1 ) 

650 FORMAT!// 10X. 'CORELLATION POINT'/ 

1 40X. 3WX.F12.3.4X) /) 

660 FORMATI///50X. 'ISOLATOR CHARACTERISTICS'/ 


-•/// 


SIDE 


= '.F12.3. 


''F12.3 / 
'.F12.3 / 

= '.F12.3 / ) 

= '.F12.3 / 

= '.F12.3 /) 


'.F12.3. 

'.F12.3 


/// 


50X. 

10X. 'PRECHARGE > HIGH SIDE = '.F12.3. 

i4x » ' air volume reaction surface > high side 

lOX.'PRECHARGE . LOW SIDE = '.F12.3. 

14X.'AIR volume REACTION surface . LOW 
lOX.'GAIN = '.F12.3. 

I4X. 'DAMPING 
670 FORMAT! 

1 10X. 'SHAKER FORCE 

2 14X. 'LONGITUDINAL ANGLE 

3 10X. 'FREQUENCY = '.F12.3. 

4 14X. 'LATERAL ANGLE 

680 FORMAT1///50X. 'PROCESSED DATA '/50X. 14 1 ' -' ) // 

1 5X» 'EXTERNAL APPLIED LOADS AT A •// 

22 X.'V • .F12.3.4X. 'S • .F12.3.4X. 'D 

3 4X ' ' MX • .F12.3.4X. 'MY • .F12.3.4X. 'MZ 

4 5X> 'CORRECTED ISOLATOR LOADS' // 

5 2X.'FI1 • .F12.3.4X. 'FI2 • .F12.3.4X. 'FI3 • .F12.3) 

690 FORMAT I///51X . ' OUTPUT ' /51X .6 ( •-•)// 

1 5X> 'MEASURING SYSTEM LOADS AT A'// 

2X. ' V • .F12.3.4X. 'S ' .F12.3.4X. 'D 
4X . 'MX ' .F12.3.4X. 'MY • .F12.3.4X. 'MZ 
5X. 'TOTAL APPLIED LOADS AT A •// 

2X»'V • .F12.3.4X. 'S • »F12.3'4X» 'D 

4X * ' MX ' .F12.3.4X. 'MY ' .F12.3.4X. 'MZ 

5X. 'Calculated external loads*// 

8 2 X.'V ' .F12.3.4X. 'S • . F12 . 3 . 4X . ' D 

9 4X< 'MX • .F12.3.4X' 'MY ' .F12.3.4X. 'MZ 
700 FORMAT!/// (6! 5X.F12. 3. 3X> /) ) 

710 FORMAT!/// (6I5X.F12.3.3X1/1 1 

720 F0RMAT11H1.59X. 'PROCESSED DATA • /60X . 14 < )/// 

1 10X. 'ACCELERATIONS: ' .25X. 'X' .17X. ' Y' » 17X . 'Z' . 15X . 'ROLL' . 14X » 

2 'PITCH' / 1 

730 FORMAT! 10X. 'UPPER PART OF UPPER BODY*/ 

1 25X. 'STEADY • .5I3X.F12. 3.3X1 / 

2 25X. 'VIBRATORY • .5I3X.F12. 3.3X> / 

3 25X. 'PHASE • .5(3X.F12.3.3X1 / 1 


» .F12.3. 
'.F12.3 // 


• .F12.3. 
'.F12.3 // 


• »F12.3. 
'.F12.3 1 


740 FORMAT! 10X. 'LOWER PART OF UPPER BODY'/ 


25X. 'STEADY 
25X. 'VIBRATORY 
25X, 'PHASE 


750 FORMAT! 10X. 'UPPER BODY' / 


25X. 'STEADY 
25X. 'VIBRATORY 


• »5<3X»F12. 3.3X1 / 

• .5I3X.F12.3.3X1 / 

• .5I3X.F12. 3.3X1 / 1 

•»5(3X»F12.3.3X1 / 

••513X.F12. 3.3X1 / 


IRBS7850 00951 
IRBS7860 00952 
IRBS7870 00953 
IRBS7880 00954 
1RBS7890 00955 
IRBS7900 00956 
IRBS7910 00957 
IRBS7920 00958 
IRBS7930 00959 
IRBS7940 00960 
IRBS7950 00961 
IRBS7960 00962 
IRBS7970 00963 
IR8S7980 00964 
IRBS7990 00965 
IR8S8000 00966 
IRBS8010 00967 
IRBS8020 00968 
IRBS8030 00969 
IRBS8040 00970 
IRBS8050 00971 
IR8S8060 00972 
IRBS8070 00973 
IRBS8080 00974 
IRBS8090 00975 
IRBS8100 00976 
IRBS8110 00977 
IR8S8120 00978 
IRBS8130 00979 
IR8S8140 00980 
IRBS8150 00981 
IRBS8160 00982 
IRBS8170 00983 
IRBS8180 00984 
IRBS8190 00985 
IRBS8200 00986 
IRBS8210 00987 
IRBS8220 00988 
IRBS8230 00989 
IRBS8240 00990 
IRBS8250 00991 
IRBS8260 00992 
IR8S8270 00993 
IRBS8280 00994 
IRBS8290 00995 
IRBS8300 00996 
IRBS8310 00997 
IRBS8320 00998 
IRBS8330 00999 
IRBS8340 01000 
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X 1234567 

0 0 0 0 0 

3 25X . 'PHASE ' .513X.F12.3.3X) / ) 

760 FORMAT!/// 10X. 'INERTIA FORCES; ' >24X» ' X' » 17X. • Y' . 17X* 'Z' »16X» 

1 'MX* <16X. 'MT' /) 

770 FORMAT!/// 10X, ' ISOLATOR LOADS; • .22X. 'NO. 1 • ' 31X. 'NO. 2* .31X. 'NO. 3' 
1 // 

2 25X> 'STEADY • >2<3X>F12.3»21X> >3X'F12.3 / 

3 25X, 'VIBRATORY • ,2(3X,F12.3,21X> .3X.F12.3 / 

4 25X. 'PHASE ' .213X.F12.3.21X) .3X.F12.3 / ) 

780 FORMAT ( 1H1 .59X. ' OUTPUT' /6 OX » 6 ( * — * ) /// 

1 10X. 'STEADY LOAOS AT POINT A! ' .5X. 'X' .16X. 'Y' . 16X. <Z> • 

2 15X. 'MX' »15X. 'MY* ,15X. 'MZ' /> 

790 FORMAT! 10X. 'MEASURED • .6X.6(3X»F12.3.2X) / 

1 10X. 'APPLIED • »6X.6(3X.F12.3.2X) / 

2 10X» 'ROTOR-MEASURED' »6X»6(3X'F12.3'2X) / 

3 10X, 'ROTOR-APPLIED • .6X.6(3X>F12.3.2X) / ) 

800 FORMAT!/// SOX. 'VIBRATORY ROTOR HEAD FORCES' // 

1 40X. ' VERTICAL (Z) ' ,23X. 'LATERAL <Y> • . 21X. 'LONGITUDINAL ! X) •/ 

2 30X.3! 3X. 'AMPLITUDE'. 11X. 'PHASE'. 6X) / ) 

810 FORMAT! 10X, 'APPLIED ' .6X.6(3X.FI2.3.2X) / 

1 10X, 'MEASURED • .6X.613X.FI2.3.2X) / ) 

C 

C 

STOP 

END 


8 

U1005 CD COUNT 

IRBSS350 

01001 

IRBS8360 

01002 

IRBS8370 

01003 

IRBS8380 

01004 

IRBS8390 

01005 

IRBS8400 

01006 

IRBS8410 

01007 

IR8S8420 

01008 

IRBS8430 

01009 

IRBS8440 

01010 

IRBS8450 

01011 

1RBS8460 

01012 

IR8S8470 

01013 

IRBS8480 

01014 

IRBS8490 

01015 

IRBS8500 

01016 

IRBS8510 

01017 

IRBS8520 

01018 

1RBS8530 

01019 

IRBS8540 

01020 

IRBS8550 

01021 

IRBS8560 

01022 

IRBS8570 

01023 

IRBS8580 

01024 


APPENDIX III-C 


Sample Test Cases 


ACTIVE ISOLATOR FULL SCALE SHAKE TEST 


PREPARED BY 
CHECKED BY 
DATE 


J.DEFELICE 
JAN. 3* 1973 


VIBRATORY EXCITATION— LONGITUDINAL 500 LB FORCE 30HZ (RUN 6C > 


MASS DATA 



UPPER PART OF 

lower part of 

TOTAL UPpER BODY 

LOWER BODY 

total AIRCRAFT 


LOWER BODY 

LOWER BODY 




w 

48o8.000 

3004.700 

7812.700 

26187.300 

34000.000 

X 

337.160 

335.650 

336.579 

336.350 

336.403 

Y 

.000 

• 000 

.000 

.000 

.000 

z 

257.120 

199.980 

235.144 

129.820 

154.022 

IX 

2060.000 

1522.000 

19222.750 

113336.000 

305494.310 

IY 

2372.000 

198Q.000 

20003.673 

1720534.000 

1913473.800 





geometry 




NO. 

X 

Y 

Z 

ANGLE 

isolators: 







1 

306.600 

.000 




2 

356.100 

25.200 




3 

356.100 

-25.200 




LOAD CELLS (MEASURING) J 
1 

324.000 

-20.000 

197.125 

2 

322*000 

18.300 

197.125 

3 

363.600 

20.000 

197.125 

4 

363.600 

-19.000 

197.125 


LOAD CELLS (TORQUE APP) : 
1 

335*700 

26.200 

197*125 

65*000 

2 

348.300 

-26.200 

197*125 

65.000 


11 6 




VIBRATORY LOAD LONGITUDINAL EXCITATION 500 LB FORCE 30HZ (Run 60 



Hh H I- O H C\J K) 

21 H (\1 H(\| K) ^ .H <\J K) C\jfO^-lOv£)r^COC7'-H-««H*H 

UJ a. Q. OO O O X J I J I J J J _l O O U DU O O O O O O U O 

»- ir<< jj j jHHWMioiouuiJuuuuuuuuuuuuu 

*-• Q »— >— li.lL U.lLQ.Q.D.Q.Q.lQQQ< <<<< <<<<<<<< 
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360.000 


I 

I 

I 

I 

I 

I 

I 

O I 

vO I 


I 

CL I 

« I 
I 

rsi I 
I I 
o I 
ro I 
I 

UJ I 
O I 
CL I 

O I 

U. I 
I 

CO I 
-I I 


I 

Z I 

O I 

M I 

h I 

< I 
I- I 
*-• I 
O I 
X I 
UJ I 

I 

_J I 

< I 

Z I 
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PROCESSED DATA 


N<00 


O 


O MO 

or*- o 

0 r* o 

•' • • 

1 i 


it- r- o 
moo 
nno 


O^O 



lor- o 
« * • 
^ lO -4 
I CM 


IO 



« 


_J 

0(0 0 
o ro o 
OHO 
• • • 

0(0 0 
o O' o 
OHO 

• • • 


OlOO 

o co o 
o O' o 

• • • 

0(0 0 
or- o 
o cm o 

_» 

o 

O 

X 

ro o 

«H O 

o 

CO 

CO 

z 

CO vO 

O v£ 

cl 




CNirO 

ro ro 



o o o 

o o o 

o o o 

*H 4 - O 


cm r* o 


ro cm o 

ro CM o 

if AO O 

if O' o 


CM tO O 


o o o 

o o o 

CM HO 

*-• O O 


HHO 


• • • 

• • • 

• • • 

• • • 


• • • 

N 

1 o 

1 O 

N f «00 

o o o 

CM 

ro in o 


AO 



O' *0 it 

• 

r*- •-* co 


ro 

ro 

-« in 

I in 

o 

OHH 




i 


z 

CO 




O «H O 

r- o o 
o o o 


o o o 
o o o 
o o o 


O O O O -H O 

vO it O OlO O 

l r> vO o O O O 



>- 


vO rO O 

ro co 

K) «H 


o 

vO 

rO 



o .f o 

o o 

o 

o v£> o 

Ot O 


>D -t o 


•H O' O 

O CM 

o 

CO o O 

ro O' o 


O' it O 


o o o 

O O 

o 

OHO 

o co o 


co r- o 


• • • 

• • 


• • • 


H 

• • • 

X 

1 1 


1 

X CO ro o 

O' o 

• 

vO CM 1 





it in co 

in co 

o 

f O 





1 it «H 

H 

z 

O' -• 
O' 



>- >- 

>- >- 


>- >- 

>- >- 


>- 


o cl 

Q CL 


O CL 

a a 


CL 


o o 

O O 


O O 

o o 


O 


cd >- r- 

CD >- »- 


cd >- r- 

CD >- K 


>- I- 


Q < UJ 

Q < UJ 


O < UJ 

Q < UJ 


Q < UJ 


cl < a: in 

cl < cl in 


CL < (L in 

QC < QC in 


< cr in 


U UJ CD < 

UJ UJ CD < 


UJ UJ CD < 

UJ UJ CD <f 


Ul CD <t 


1 »-hX 

o. r— »-* i 


Q. h-Ml 

Q.I-HI 


»- X 


a. in > o. 

a in > a. 


a. in > o. 

CL in > CL 


in > o_ 


Z> 

3 

•• 

z> 

3 

•• 


• • 



in 



in 


in 

u. 

IL. 

!J 

U_ 

u. 

a 


z 

o 

O 

o 

o 

o 

< 


o 



CL 



o 


M 

t- 

»- 

O 

»- 

r- 



r- 

cl 

CL 


CL 

CL 



< 

< 

< 


< 

< 

(L 


a: 

CL 

CL 

< 

0. 

CL 

O 


lU 



*— t 



»- 


_J 

cl 

CL 

h 

CL 

CL 

< 


u 

UJ 

UJ 

CL 

UJ 

UJ 

_l 


o 

CL 


UJ 

CL 

* 

o 


o 

Q. 

O 

Z 

CL 

O 

in 


< 

3 

_l 

H 

3 

-1 

»-« 
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H UJ • • 

(DO o sO 

Z 3 O (M 

o i- in m 


CM * * CM CM -H ^ 


O <0. 

_J O LlI Q_ 

UJ QI< 

>- tr uj i i 

O 3 ~0C CE 

< l/l JOO 

UJ < Q_»— I— 
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APPENDIX IV 


Analysis of Contributions to Error in 


Measuring Vibratory Loads 


The following is an analysis of three vibratory load conditions which sub- 
stantiate the major contribution of acceleration measurements to the determina- 
tion of vibratory loads . 

1 . Vertical Vibratory Load - 500 lbs. @ 30 Hz. 


Isolator Vibratory Loads 
#1 10U.5 lbs. 

#2 90.5 Its. 

ft 3 64.0 lbs. 


Total 259.0 lbs. 

Inertia Response of Transmission/Rotor Head 

Acceleration - g's Force 


Rotor Head 
Transmission 


.06 


288 lbs. 

.06 


180 lbs. 


Total 

468 lbs. 


% 


Total measured load = 727 lbs. 

Accelerometer contribution to Total measured load = 468 lbs. 


2 . Vertical Vibratory Load - 495 lbs. % 30 Hz . 
Isolator Vibratory Loads 


# 1 -102 lbs. 

#2 +102 lbs. 

#3 -116 lbs. 

Total -11 6 lbs . 

Inertia Response of Transmission/Rotor Head 
Acceleration - g's 
Rotor Head .06 

Transmission .06 

Total 


Force 


288 lbs. 

l80 lbs. 
468 lbs. 
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Total measured load = 352 lbs . 

Accelerometer contribution to Total measured load = 468 lbs . 


Longitudinal Vibratory Load - 1000 lbs. @ 18.5 Hz 


Loads present in Longitudinal Drag Links 

Forward +15 Its. 

Aft - 6 5 lbs . 

Total -50 lbs . 

Inertia Response of Transmission/Rotor Head 
Acceleration - g's 
Rotor Head .23 

Transmission .07 

Total 

Total measured load = 1270 lbs. 


Force 

1105 lbs. 

210 lbs. 
1315 lbs. 


Accelerometer contribution to Total measured load = 1315 
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